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I. SUMMARY AND INTRODUCTION 


1 . 1 SUMMARY 


The VTOL Approach and Landing Technology (VALT) program is a 
significant experimental research program aimed at establishing a 
data base for rotorcraft operation in a terminal area environment. 
The work described in this report was undertaken to determine heli- 
copter math models suitable for analyzing maneuvers along a VTOL 
trajectory and to apply these math models to determine the flight 
test procedures of greatest effectiveness in establishing helicop- 
ter dynamic characteristics in this mode of operation. As the 
principal result of this investigation, a flight test specifica- 
tion is presented for the CH-47 VALT aircraft operating along the 
specified VTOL trajectory of the VALT program. 

Simulation mathematical models are an integral part of the 
process that led to the recommended flight test specification. 

These models are used to evaluate possible error sources and to 
validate input design recommendations. Two CH-47 simulations were 
used for this study: a simulation based on CH-47 linear derivative 

models reported in Ref. 1 and the C-81 Rotorcraft Flight Simula- 
tion Program Ref. 5 (AGAJARMY Version). The instrumentation 
error models and LRC data recording system error contributions 
were simulated as in Ref. 3. 

The principal results of this study are summarized as 
follows : 


(I) Specification of flight tost procedure for parameter 
identification of the CH-47 VALT Research Aircraft 
along a VTOL trajectory representative of VALT pro- 
gram requirements. This specification includes 
requirements for: 


• Test points 



• Control input shapes 

• Stability augmentation system gains 

• Maximum wind gust intensities for acceptable 
accuracy 

(2) Investigation and determination of the influence on 

derivative identification accuracy of: 

• The CH-47 force-moment characteristics along a VTOL 
trajectory, linear model representation, rotor degrees 
of freedom and SAS characteristics 

• The CH-47 VALT Research Aircraft instrumentation and 
the LRC data recording system 

• Control input design and procedures 
1.2 INTRODUCTION 

1.2.1 VALT Program Overview 

The NASA Langley Research Center (LRC) VTOL Approach and 
Landing Technology (VALT) Program encompasses requirements 
analysis, design studies, systems development, ground simulation 
and flight validation of VTOL aircraft, subsystems and tech- 
nology developments. The VALT research aircraft, a CH-47 tandem- 
rotor medium- transport helicopter, is used to provide the neces- 
sary system design evaluations and flight test validation of 
fully integrated system developments. A photograph of the 
LRC VALT CH-47 Research Aircraft is shown in Figure 1.1. 

One VALT system development of significance is associated 
with the control system design for approach and landing condi- 
tions. Such control design is dependent on a model of the 
vehicle which can be used for specification of gains to improve 
vehicle handling qualities, stability, or overall integration with 
the guidance systems. To achieve such a model, the technology of 
system identification is to be applied to VALT CH-47 flight test 
data . 




Figure 1.1 The LRC CH-47 VALT Research Aircraft 



1.2.2 Study Overview 


Purpose . The purpose of this investigation is to produce 
a flight test specification for the LRC CH-47 VALT Research 
Aircraft for parameter identification along the VTOL trajectory 
of the VALT program. In meeting this objective, specific areas 
of investigation are addressed in this report, which include: 

(1) Selection of a nonlinear simulation program of 
the CH-47 and sufficient test points along the 
specified VTOL trajectory; 

(2) Selection of an appropriate linear model represent- 
ation of the CH-47 along the VTOL trajectory (e.g. 
test points requirements) ; 

(3) Analysis of control input designs for parameter 
identification, and consideration as to the feas- 
ibility of using the CH-47 production stability 
augmentation system. 

Scope and Limitations . The results presented in this report 
represent a rotorcraft parameter identification simulation study, 
which includes a systematic investigation of sources contributing 
to error in parameter identification. The results include the 
instrumentation error methods of Ref. 3. These results contribute 
to the VALT program objectives by providing the necessary analysis 
leading to a flight test specification for accurate CH-47 param- 
eter identification. The CH-47 flight test specification result- 
ing from this investigation is specifically designed for extracting 
the stability and control coefficients necessary for the VALT 
program control system design objectives. 

Parameter identification accuracy is related to the quality 
of the control input design specification, which is dependent on 
the fidelity of t he m athematical model used for the synthesis. 
Several factors influence the mathematical model selection: 

(1) Degrees of freedom 
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(2) Significance of coupling between degrees of 
freedom 

(3) Adequacy of a linear vehicle representation (and 
a nonlinear model representation). 

These considerations are particularly important for rotorcraft 
because of their characteristic high order models (rotor and 
fuselage) , inherent coupling between fuselage longitudinal and 
lateral modes in some vehicles, and the uncertain effects of 
aerodynamic nonlinearities under some flight conditions. The 
attention which must be given to each of these factors depends 
most importantly on the particular rotorcraft, its rotor type 
and its stability augmentation system. For the VALT aircraft 
and the objectives of the VALT study, a linear model is speci- 
fied. The extent of this effort includes the following: 

(1) Determine an adequate linear model approximation 
to the CH-47 vehicle equations of motion (based 
on nonlinear simulation evaluation) . 

(2) Design inputs for these linear models (as well as 
specific test points to be tested) . 

(3) Evaluate the anticipated accuracy obtainable from 
the designed inputs and specify the post-flight 
data processing requirements necessary to account 
for errors in the assumed input design model. 


Approach . The technical approach is graphically outlined 
in Figure V.2, and a summary of each element of the figure is 
discussed below. The following discussion summarizes the 
technical approach and the rationale leading to the procedure 
followed. 


The CH-47 nonlinear simulation used in this analysis (indi- 
cated in Figure 1.2) is the AG A J ARMY version (1975) of the C-81 
Rotorcraft Flight Simulation (Ref. 5) computer program developed 
by Bell Helicopter Co., Fort Worth, Texas, under direction of the 
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Figure 1.2 Technical Approach 




Eustis Directorate, USAAMRDL, Fort Eustis, Virginia. The C-81 
computer program is operational on an IBM 360/65 computer from 
which all of the C-81 results presented in this report were 
obtained . 

The C-81 program is used in two modes of operation in this 
study: (1) in a nonlinear time history simulation mode, and (2) in 

a trim and stability derivative mode. The original plan for 
this investigation was to use the nonlinear C-81 time history 
simulation mode with parameter identification which would sub- 
sequently validate the linear simulation identification. The 
computer run time associated with a single nonlinear C-81 time 
history simulation was found to be prohibitive for the data 
lengths required for parameter identification. For example, simu- 
lating 3.5 seconds of real time data required 30 minutes of CPU 
t ime . 

, » 

Since data record length requirements for identification are 
typically 15, seconds, this would require in excess of two hours 
of computer CPU time for a single maneuver. Thus, it was decided 
to use short time history data lengths of 0.5 second from the 
nonlinear simulation model for comparison purposes to the various 
linear models used. The linear 10 degree-of -freedom (DOF) rotor/ 
fuselage model (obtained from the stability derivative option 
of C-81) was subsequently found to provide an excellent time his- 
tory comparison to the nonlinear C-81 and, thus, the 10 DOF 
linear model was used as the baseline model throughout this 
inves tigation . The C-81 stability derivative option provides 
for both a 10 DOF derivative model with separate rotor flapping 
degrees of freedom and the conventional 6 DOF quasi-static 
derivative array. Both models (10 DOF and 6 DOF) are used in 
this study. 
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The linear models developed from the C-81 simulation 
program consist of: (1) the 10 DOF fuselage/rotor model with 
first harmonic rotor representation, (2) the 6 DOF quasi- 
static derivative model, (3) a 10 DOF fuselage/rotor model with 
first-order rotor representation (flap advancing mode deleted) , 
and (4) the above models with the 8-state variable CH-47 produc- 
tion stability augmentation system (SAS) . The C-81 stability 
derivative option provides derivatives in the form of a second 
order model (mass, spring, damper) and requires conversion to 
state variable form. The three aircraft models in state variable 
format are obtained from the computer program QUASI -STATIC , which 
accepts as input the C-81 second order derivative equations and 
provides the three models in state variable form as output, along 
with eigenvalues and eigenvectors for each mo del . 

The derivative models, taken from Ref. 1, were used mainly 
for comparisons with the C-81 values. A detailed comparison 
was made and it was found that the speed-dependent derivatives 
(e.g. M u , M w , etc.) showed large discrepancies between the 
Ref. 1 and C-81 derivatives, and yet the damping and control 
derivatives showed adequate agreement. The speed-dependent 
derivative discrepancy was traced to the fact that the C-81 
model does not contain rotor- to-rotor interference effects; 
whereas, the model taken from Ref. 1 does. Because of the 
importance to handling qualities and control system designs of 
the speed-dependent derivatives, it was decided that the deriv- 
ative model taken from Ref. 2 is more representative of the 
CH-47 and should be used as the primary simulation, whereas 
the C-81 models should be used primarily to investigate the 
effects of rotor degrees of freedom on derivative identification 
accuracy. Thus, both simulation models are used for this study. 

The control input designs indicated in Figure 1.2 are based 
on the integrated application of eigensystem analysis (e.g. 
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location of system time constants) and iterative evaluation of 
simulated responses of the vehicle. 

The data generation and wind gust simulation indicated 
in Figure 1.2 use the appropriate linear model (either C-81 derived 
or Ref. 1 model) and the designed control inputs. A linear simula- 
tion computer program is used for this purpose. The generated data 
are stored on magnetic tape and provide the data base for further 
analysis. Wind gust is simulated using a simplified Dryden model. 

The CH-47 instrumentation and LRC data acquisition system are 
modeled as in Ref. 3. The magnitudes of the error sources are ob- 
tained from the existing CH-47 instrumentation specification 
sheets, the LRC-piloted aircraft data system (PADS), and the LRC 
Research Aircraft Ground Station Data Reduction System (RAGS) spec- 
ifications. ; 

The parameter identification algorithm used in this investi- 
gation is the Systems Control, Inc. Maximum Likelihood Identif i- 
cation program (SCIDNT) . This program was. used in identification 
of open-loop stability and control coefficients with the CH-47 
production SAS included as part of the model. In addition, the 
Optimal Subset Regression (OSR) is used for model structure 
determination. 

The flight-test specification for the CH-47 VALT Research 
Aircraft is obtained after simulations using the analysis methods 
summarized above and indicated in Figure 1.2. Flight test points 
presented in Chapter IV are based upon both quantitative and 
qualitative results. 

1.2.3 Report Organization 

Subsequent chapters of this report are organized as follows. 
Chapter II describes the factors which influence parameter identi- 
fication accuracy, with specific emphasis given to the CH-47 



characteristics along the VALT trajectory. Chapter III presents 
the simulation results which formed the basis for the CH-47 flight 
test specification for parameter identification. The recommended 
parameter identification flight test plan is presented in Chapter 
IV. The conclusions and recommendations are given in Chapter V. 



II. PARAMETER IDENTIFICATION ERROR SOURCE ANALYSIS 


This chapter describes those factors which influence para- 
meter identification accuracy, with specific emphasis on the 
CH-47 along a VALT trajectory. The physical considerations which 
are of fundamental importance for accurate identification are dis- 
cussed. It is often assumed that an accurate parameter identifica- 
tion algorithm is all that is required for accurate identification. 
The accuracy in any parameter identification attempt is, however, 
closely related to knowledge of the physical characteristics of the 
system to be identified, the environment, the instrumentation 
measurement system, and the control input including constraints. 
These factors are discussed since they are of paramount importance 
to the task of defining a flight test specification. Each section 
in this chapter addresses in detail a major factor influencing 
parameter identification accuracy. These include: (1) the VALT 

trajectory, (2) the CH-47 characteristics and SAS, C3) the CH-47 
linear model representation, (4) instrumentation and LRC data 
recording facilities, (5) the flight test environment, and (6) 
flight test inputs. Following this, the next chapter addresses 
the important considerations in the selection and use of parameter 
identification algorithms and presents pertinent theoretical 
developments . 

2.1 THE VALT TRAJECTORY 

i 

The flight test specification resulting from this investi- J 
gation for parameter identification is to be along a nominal j 
VALT trajectory. A sufficient number of flight test points are | 
required to permit subsequent validation and possible improve- 
ment of CH-47 simulations along this trajectory. The stability 
and control characteristics of helicopters in general and. 
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specifically , the CH-47, vary greatly over the flight envelope 
and are particularly dependent on velocity (longitudinal, verti- 
cal, and lateral). Thus, the particular nominal flight path along - 
which testing is to occur will result in considerable variation 
of the stability and control derivatives and will also determine 
the degree of aerodynamic and dynamic nonlinearity at each test 
point. In addition, the proximity to the ground (in-ground 
effect) and steady wind shear profile penetration may be relevant 
depending on the flight paths. 

Figures 2.1 and 2.2 define the nominal VALT approach trajectory 
which was used for this study. Important transition points are 
indicated by the letters A through D. The curved flight path from 
points A to C is flown at a constant airspeed (65 knots) . Altitude 
is held constant at 229m ft between points A and B. At point B, 
a 2.54 m/sec descent rate is initiated and maintained through 
point D. At point C, a velocity deceleration to hover is initi- 
ated while altitude decreases from 137m to 15m. 

The shape of this flight profile is dictated by the many 
constraints imposed on commercial transport VTOL operations in 
the terminal area, which are discussed in general in Ref. 9. The. 
area of main emphasis for parameter identification is between 
points C and D--the terminal descent approach with decelerating 
velocity. This transition region, which is schematically illus- 
trated in Figure 2.3, is the most difficult to model with analytic 
simulations, usually contributes the most adverse handling quali- 
ties, is the region of least stability, and is thus most demanding 
on control system designs and navigation in minimizing flight 
path dispersions. 

2.2 CH-47 CHARACTERISTICS 

The CH-47 is a twin-engine, tandem-rotor helicopter designed 
for all-weather, medium-sized transport operations. The rotors 
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Figure 2.1 VALT Mission Ground Track 
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Figure 2.3 VALT Transition Flight Profile 



each contain three fully articulated blades with pitch, flap, and 
drag hinges. The rotors are driven in opposite directions and 
primary helicopter control is obtained as follows: (1) longitudinal 

control (pitch) by differential collective, (2) vertical control 
by ganged collective, (3) lateral control (roll) by later cyclic, 
and (4) directional control (yaw) by differential lateral cyclic. 
This section discusses those characteristics particular to the 
CH-47 which are of importance to parameter identification. A gen- 
eral discussion of the CH-47 (Models A, B, and C) is presented in 
Ref. 12 with emphasis on general stability and control. Reference 
13 discusses the ground-based flight simulation of the CH-47C, and 
Ref. 1 presents the equations and model to be used to simulate the 
CH-47 VALT Research Aircraft which will be flight tested for param- 
eter identification. 

2.2.1 SAS-Off Stability and Control Characteristics of the 
CH-47 


This section discusses the basic open-loop stability and 
control characteristics of relevance to parameter identification. 
These characteristics are obtained from the various analytic 
models presented in Chapter III and Appendices A and B, and 
are introduced here to emphasize those factors influencing the 
flight test specification. The basic trend in stability and 
control characteristics of helicopters, in general, is similar 
for low-speed flight and the significant differences in handling 
qualities are a result of the magnitude of these characteristics. 


Figure 2.4 presents the CH-47 characteristic root locations 
with SAS-off for two airspeeds (Vertical descent and 40 kts) at 

? C * -n ' a .i „ _ „ - — -- ^ ^ - v ■* — 1 3 c i — t r> /-> i • 3 r>r.i a n j „ 

i. a 1 * ul/ act ucatcilL lace lUuiaincu x. i out cue l/Ao-sujjpiiCvi / ae “ 

rivative model). At hover, both the Dutch roll and longitudinal 
long-period modes are shown to be slightly unstable and oscillatory 
(refer to Fig. 2.4, open triangle symbol). As airspeed increases 
to 40 knots, the Dutch roll remains nearly fixed, whereas the 
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Figure 2.4 CH-47 SAS-Off Characteristic Root Location with Increasing Speed 




longitudinal long-period mode migrates further toward the right 
half plane and is unstable and aperiodic at a value of +0.5. 

Figure 2.4 shows that the stability characteristics change con- 
siderably from hover to 40 knots. This variation in vehicle 
dynamics results from the influence that airspeed has on the 
stability parameters, which make up the equations of motion. 

Figure 2.5 illustrates the type of nonlinear characteristics 

which must be considered for rotorcraft modeling beyond the 

linear regime. This figure presents predicted speed-dependent 

CH-47 derivatives M and M as a function of airspeed for 

u w 

three values of descent rate (h = -5.08, -2.54, 0 m/sec). The 
plots of M and M presented in this figure emphasize three 
main considerations. First , these derivatives are very nonlinear 
with airspeed and thus the perturbations in airspeed must be 
kept typically less than +_5 knots. Second , these derivatives 
also are very nonlinear with descent rate and thus perturbations 
in la must also be kept small. Third , at hover and 20 knots 
airspeeds, in order to retain approximate linearity in M u , 
vertical velocity must be maintained near its trim value (i.e. 
negligible or near-zero perturbation). Likewise at hover, in 
order to retain linearity in M , airspeed must not be allowed 
to deviate by any substantial amount from its trim condition (a 
+5 knot change in airspeed causes to change by 0.005 at 

hover) . 

The above discussion indicates that in order to identify 
the approximate linear stability derivative M y at hover and 
20 knots, perturbations in li should be kept to a minimum and 
only airspeed should be perturbed. In an analogous manner, the 
identification of M must be done by perturbing h while 
maintaining nearly constant airspeed. Thus, the design of control 
inputs must be such that they contain state variable constraints. 
The reader is directed to Appendix D for a discussion that valid- 
ates the qualitative conclusions of this paragraph. 
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Figure 2.S CH-47 Nonlinear Stability Characteristics at Low Speed (data from 

ref. [1]) 


2.2.2 SAS-On Dynamic Characteristics of the CH-47 


The previous section reviewed those characteristics of the 
basic vehicle (SAS-off) which are of concern to parameter identi- 
fication. Because the vehicle is unstable without a SAS , the 
CH-47 production SAS is investigated in this study as to its 
effect on parameter identification. 

Figure 2.6 shows the CH-47 characteristic root location with 
SAS-on at hover and 40 knots airspeed and 2.54 m/sec descent rate 
At the hover condition, both the Dutch roll and longitudinal long- 
period modes are stabilized by the SAS (open triangle symbols) , 
whereas at the 40 knot condition, only the Dutch roll is stable 
The longitudinal long-period root i- still unstable, although with 
a greater time to double amplitude. The frequency of the pitch 
and roll short-period modes is found to be faster with SAS (as 
seen by a comparison of the pitch and roll short period roots 
presented in Figures 2.4 and 2.6). 

Before completing the discussion of SAS-on dynamic character- 
istics, it is necessary to briefly review the effect that the rotor 
degrees-of-freedom have on vehicle dynamics. Figure 2.7 shows the 
approximate location of basic rigid body characteristic roots 
both with SAS-on and SAS-off. Rotor characteristic root location 
for three different helicopters are on this figure to emphasize 
the proximity of the CH-47 rotor mode in comparison to the CH-53 
and XH-59A (ABC) helicopters. The rotor root (flap regressing 
mode) of the CH-47 is shown to be far removed from rigid body 
modes. The results reported in Ref. 2 indicated the importance 
of the rotor degrees of freedom on identification of rigid body 
derivatives and presented results for the CH-53. The frequency 
of the CH-47 rotor flap regressing mode shown in Figure 2.7 is 
greater than that of the CH-53 and thus, the influence of the 
rotor should be correspondingly less. It should be noted that 
the rotor roots shown in Figure 2.7 are obtained with the SAS-off 
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Figure 2.6 CH-47 SAS-On Characteristic Root Location with Increasing Speed 
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Figure 2.7 CH-47 Rotor Roots Are Far Removed from Rigid Body Modes 




and, thus, ignores SAS-rotor — fuselage coupling. The CH-47 pro- 
duction SAS does couple with the rotor degrees -of -freedom . This 
coupling is shown in Figure 2.8, which shows the rotor flap re- 
gressing mode with SAS-off, CH-47 production SAS, and CH-47 pro- 
duction SAS with roll rate feedback gain divided by two. 

A reduction of the SAS roll rate feedback gain from its 
production value of = -12.7 to = -6.35 cm/rad/sec reduces 
the rotor fuselage coupling as indicated by the triangle symbol 
of Figure 2.8. To minimize rotor-rigid body coupling, the CH- 
47 flight test identification should be conducted with the 
CH-47 production SAS operative with the roll rate feedback gain 
halved. The identification results presented in the chapters 
to follow contain this SAS modification. 

The use of stability augmentation allows independent excit- 
ation of vertical velocity with collective input and independent 
excitation of longitudinal velocity with longitudinal stick input. 
The ability to perform this independent velocity excitation is 
a direct result of the SAS which, by nature of the control law 
structure, produces pole-zero cancellation. The high degree of 
decoupling will be shown by time history responses in a later 
section. This decoupling can be an advantageous feature for 
input design. 

2.3 MATHEMATICAL MODEL STRUCTURE 

This section discusses the development of the linear models 
from C-81 and Ref. 1 models, the linear model structure used for 
the simulation, the linear model structure used in the identifi- 
cation and the CH-47 production SAS mathematical representation. 

As discussed previously, the models developed from C-81 and those 
supplied by Ref. 1 are both used in this analysis. The analytical 
form and coefficient values of these models are presented in 
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Figure 2.8 Effect of Roll Rate Feedback on Rotor -Airframe Coupling 




Appendices A and B. Before discussing details of the linear 
models, a discussion of the results and developments from C-81 
is presented which includes time history comparisons of the 
linear models and the nonlinear C-81 simulation. 

2.3.1 Nonlinear C-81 Simulation, Linear Model Developments, and 
Time History Comparisons 

The nonlinear C-81 simulation program is used to define trim- 
med conditions from which linear perturbation models can be 
determined or from which time histories can be generated. Both 
of these modes of operation are used and discussed below. 

Models Derived from the Nonlinear C-81 Program . Two models 
are derived from the nonlinear C-81 computer program: viz, a 10 
DOF model and a 6 DOF quasi -static model. The 10 DOF model is 
obtained from C-81 by perturbing each independent variable sep- 
arately while keeping all others at their value, computing the 
vehicle force and moments, and dividing by the independent vari- 
able. The form of the output is that of a second-order matrix 
equation as 

MX + DX + KX = BU (2.1) 

where 

M = mass matrix (10 x 10) 

D = damping matrix (10 x 10) 

K = stiffness matrix (10 x 10) 

B = control matrix (10 x 4) 

The independent variables are represented by the vectors X, X, 
and X, and the control vector by U. These are defined 
as 
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The approximate 6 DOF quasi-static derivative model is obtained 
by perturbing only the fuselage independent variables and allow- 
ing the rotor degrees of freedom to reach steady state. The 
effect of the rotor is thus lumped into the fuselage resulting 
in quasi-static derivatives. The form of the derivative equation 
is the same as shown in Eq. (2.1), except that the matrices are now 
of dimension 6x6. The independent variables are defined 
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Models Derived from the 10 DOF Linear Model . Another model 
is derived from the 10 DOF linear model-ra 6 DOF quasi-static 
model. This model is developed from the 10 DOF model by first 
converting this model to state variable form. The state variable 
representation for the 10 DOF model is 


x 


Fx + Gu 


(2.4) 



where 


F = 17x17 matrix of stability derivatives 
G = 17 x 14 matrix of control derivatives 

and the state vector x is defined in Eq. (2.5) 

x T = [e^Yuwqvpr ai MR bi MR a lTR (2.5) 

t>l T R e 1mr b lMR a lTR ni TR^ 

T 

u = ^ 6 coll 6 long 6 lat 6 rud^ 

The state model Eq. (2.4) can be rewritten and partitioned for 
development of an approximate quasi-static model. The parti- 
tioned form is 


where 



( 2 . 6 ) 


(2.7) 


The quasi-static condition occurs when the rotor state variables 
reach steady state (e.g. X 2 and X 2 are zero). Making these 
substitutions into Eq. (2.6) and solving for X 2 results in 


X 2 


F 23 *- F 21 X 1 + G 2 U ) 
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( 2 . 8 ) 



Substitution of Eq . (8) into the equation of Eq. Q2.6) results 

in the final form for the quasi-static model as 


= CF 11 - F 13 F 23 F 21> X 1 + ( G 1 ‘ F 13 F 23 G 2^ u C2 ' 9) 


QS 


3 QS 


Equation (2.9) defines the 6 DOF quasi-static model derived from 
the 10 DOF linear model. This model can be compared to the quasi- 
static model obtained directly from the C-81 nonlinear computer 
program. The state variable form for the C-81 6 DOF quasi- 
static model is 


x ' F C-81 X + G C- 81 U 


<€= 6 DOF C-81 Model (2.10) 

obtained by per- 
turbation 


where the notations F c _ gl an< * ^C-81 indicate that this model 
was obtained directly from the nonlinear C-81 program by perturb- 
ation of the nonlinear model. The derivative values for F^g-^ 
are found to be of nearly the same value as F^g obtained from 
Eq. (2.9). However, the control derivatives do not agree 

with the values computed by Eq. (2.9). This is because Gpg com- 
puted from Eq. (2.9) is derived from the 10 DOF model, which was 
found to contain incorrect control derivatives and is in error due 
to a spike anomaly occurring by perturbing at t = 0 in the C-81 
nonlinear computer program. 

The anomalous behavior associated with the C-81 computer 
simulation has been observed previously. Whenever sharp edge pilot 
control inputs were used, spike responses occurred just as when per- 
turbing the nonlinear model to derive the 10 DOF control derivatives. 
Analysis of this problem indicates that the spike results from the 
generation of the fuselage control derivatives of the 10 DOF 
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derivative model. The 10 DOF model is derived by perturbing the 
nonlinear model at t = 0 . When this is done, spike responses 
result at t = 0 and, as a consequence, the t = 0 derived control 
derivatives are in error. This error can be corrected by using 
the 6 DOF C-81 quasi-static derivatives and reversing the procedure 
employed to develop quasi-static derivatives. The 10 DOF model 
used in this study and referenced in the remainder of this report 
is the corrected 10 DOF model. 

The spike anomaly found in the C-81 program has been ob- 
served in other rotorcraft simulation programs. For example, the 
Sikorsky Aircraft General Helicopter Simulation program has been 
found to contain spike anomalies of a similar nature. Reference 

15 (page 184, Figure 1) shows the stability derivative Zt, 

a lS 

obtained by perturbing the nonlinear simulation. At time t=0 
in this figure, a large spike input exists which is an incorrect 
representation of the 9 DOF derivative value discussed in this 
reference. 

The errors in the derived control derivatives G^g resulting 
from the spike anomaly can be corrected since it is known that 
the 6 DOF quasi-static derivatives, G^g^, obtained from C-81 
are correct (since they are not obtained at t = 0, but after 
rotor transients subside). This correction follows since the 
control derivatives of Eq. G.2 . 9 ) can be set to the C-81 values as 
shown in Eq. (2.12). 


G 1 ‘ F 13 F 23 G 2 " G C- 81 


( 2 . 11 ) 


Therefore, the -10 DOF control derivative values, G-^ , are 
modified as shown in Eq. (12). 


G 1 = G C- 81 


F 13 F 23 G 2 
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( 2 . 12 ) 



where G^ represents the modified correction. This corrected 
control derivative matrix can then be used in place of 
in Eq. (2.6) to yield a 10 DOF model which does not contain errors 
resulting from the spike anomaly from C-81. This model will then 
reduce to the correct 6 DOF quasi-static derivative model of 
Eq. (2.10). 

The corrected 10 DOF linear model (i.e. with Gp is now 
a consistent representation to the nonlinear C-81 simulation 
program. This correction is justified, since the 10 DOF must 
reduce to the C-81 quasi -static model. The correction thus 
eliminates the errors resulting from the erroneous spike and 
guarantees reduction to the correct quasi-static derivatives. 

The time histories to be shown subsequently do demonstrate 
this validity. 

The above discussion presented the development of the 6 DOF 
quasi-static model from the 10 DOF (C-81 derived) linear model. 

The control derivatives of the 10 DOF were found to be in error 
and corrected in the manner discussed above. 

It is of interest to compare the eigenvalues for the two 
models--10 DOF first harmonic rotor and 6 DOF quasi - static . 

Table 2.1 shows this comparison for the hover condition and 
excellent agreement exists for the low frequency rigid body 
modes. (If this agreement did not exist, then the quasi-static 
assumption would not be a valid approach for low frequency 
modeling) . 


C-81 Time History Comparisons (Nonlinear vs. Linear) . Time 
history comparisons were conducted using the nonlinear C-81 com- 
puter program and results of the various linear models. As dis- 
cussed previously, the computer run times associated with the 
nonlinear C-81 time history simulation precluded its use for 
derivative identification. The approach taken is to validate 
that the linear simulation models adequately match the 
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Table 2.1 

Eigenvalues from the Two Linear Models (Hover) 


MODE 

C-81 
10 DOF 

REDUCED 
6 DOF 

Long. 

Divergence 

.2921 

.2982 

Dutch Roll 

.0696+j . 3920 

. 0768+j . 3850 

Vert. /Long. 

-. 1447+ j . 2935 

-.1422+j. 2973 

Pitch 

Damping 

-.9377 

-.9217 

Roll 

Damping 

-.8305 

-.8481 

Spiral 

-.4015 

-.3967 

Flap 

Regressing 

- 18 . 24+ j 6 . 381 
-17.83+j7.691 


Flap 

Advancing 

- 18. 20+ j 44. 99 
-18.31+J43.92 



' % ~~ r -7 r 


nonlinear simulation, thus justifying their use for generating the 
data for use in the derivative identification. 

Figure 2.9 presents the time history responses of the CH-47 
at the 40 knots condition due to a 1.27 cm by 0.3 sec longitudinal 
pulse input (differential collective). Figure 2.9 shows pitch 
acceleration generated from the nonlinear C-81 simulation program. 
The nonlinear simulation was terminated after 0.5 sec due to the 
large execution time associated with the C-81 program. This 
amount of simulation time is sufficient for time history comparison 
since most of the transient discrepancies occur in the vicinity 
of the pulse input. The spike anomaly discussed in detail 
previously is shown on this figure. The simulation results of 
the linear 6 DOF quasi-static model and the 10 DOF fuselage 
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Figure 2.9 Longitudinal Time History Comparison of the 
Nonlinear C-81, 10 DOF, and 6 DOF, Model 
Response (v = 40 knots, S T ~ xir , = 0.5° for 0.3 
seconds) 
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rotor model are also shown. The 10 DOF model is shown to 
represent the transients of the nonlinear response with the 
exception of the erroneous spike. The 6 DOF response does not 
model these transients; however, after the input is removed, 
the 10 DOF and 6 DOF show excellent comparison. The comparison 
between the 6 DOF and 10 DOF models is considered quite good in 
the vicinity of the pulse input. This is because pitch acceler- 
ation is primarily a result of differential collective and thus 
does not have a large time delay as would be found if pitch 
was produced by rotor flapping. 

Figure 2.10 shows the time history response due to a 1.27 cm 
in by 0.3 sec lateral pulse input (lateral cyclic). Roll 
acceleration from the nonlinear C-81 simulation, the linear 6 
DOF quasi-static, and the 10 DOF fuselage-rotor response are 
shown. A 0.2 sec delay is seen in Figure 2.10 due to rotor 
flapping. This delay is also shown in the 10 DOF linear response. 
The 6 DOF quasi-static response does not reproduce this delay • 
since the rotor degrees of freedom are not included separately 
in this model. 

Figure 2.11 shows the cross -coupling effect of the lateral 
cyclic input on pitch acceleration from the nonlinear C-81 
simulation. The response of the 10 DOF and 6 DOF linear models 
due to the lateral input is shown in Figure 2.11. The linear 
10 DOF cross-coupled response is shown to accurately represent 
the nonlinear C-81 model. The 6 DOF quasi-static response, 
however, does not represent the cross-coupling response. 

The comparisons presented in this report are represent- 
ative of the remaining time history comparisons which have 
not been shown. The comparisons shown in Figures 2.9 through 
2.11 are presented since these time histories demonstrate the 
main features of the various models which are used in this study. 



P-OOf <OEG/SEC««2J P-OQT J0EG/SEC»*2J 




Figure 2.10 Lateral Time History Comparison of the Nonlinear 
C-81, 10 DOF, and 6 DOF Model Responses (v = 40 
knots, f° r 0*3 seconds) 
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Figure 2.11 Cross-Coupled Time History Comparison of 
the Nonlinear C-81, 10 DOF, and 6 DOF 
Responses (v = 40 knots, S T = 0.5° for 
0.3 seconds) 
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The treatment of nonlinearities due to large state vari- 
able perturbations from trim was discussed in the previous 
section (SAS-off stability and control characteristics of the 
CH-47). The nonlinear trend due to speed effects obtained 
from the C-81 program is not representative of the actual vehicle 
due to the lack of inter-rotor interference modeling in this 
.program. Thus, nonlinear effects are treated based upon the 
Ref. 1 models and C-81 derived models were used primarily for 
assessing the influence of rotor degrees of freedom. 


Reference 1 CH-47 Derivative Models and Comparison to C-81. 
The CH-47 derivative models employed in this study were obtained 
from the NASA Langley Research Center and are discussed in 
detail in Ref. 1. Numerical values of the 6 DOF derivatives are 
presented in the appendix of Ref. 1 at various forward velocity 
trim conditions and descent rates. The Ref. 1 derivative models 
used in this identification effort are presented in Appendix B. 

The C-81 derivatives discussed in the previous section are 
compared with those taken from Ref. 1. This was done for three 
speeds: hover, 20 knots, and 40 knots. Table 2.2 presents a 

detailed comparison of the 40-knot condition. Those derivatives 
usually considered important are noted with an asterisk. The 
remaining derivatives are cross -coupling derivatives or contri- 
bute a negligible influence on vehicle response for small per- 
turbations and thus large differences may exist but they are of 
negligible consequence. 


The derivative correlation at hover and 20 knots airspeed is 
equally good with the exception of M u and M jt . These deriv- 


w 
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shows a considerable discrepancy to the Ref. 1 values. Table 2.3 
presents M and M for the three conditions. 


u 


w 
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Table 2.2 CH-47 Derivative Comparison (Ref. 1 and C-81) 

(V = 40 Knots, ft = -2.54 m/s) 



<U 

> 


to 

> 

0) 

o 

O CO 

OC I 
— J CJ 

I I ' 

u s- 
0) a i 

11 

c c 

TJ 
Ut C 
*— * C\J 



* 


37 












































T ab 1 e 2.3 

Comparison of and from Ref. 1 

and C-81 


AIRSPEED 

M 

u 

M 

w 

LRC 

C-81 

LRC 

C-81 

Hover 

.01356 

-.000457 

.00285 

.000104 

20 Kts 

.00965 

-.00117 

.01226 

.00162 

40 Kts 

-.00232 

-.00215 

.02073 

.00370 


The general conclusions concerning the correlation of 
derivatives from C-81 and the Ref. 1 program are that for the 
most part, good agreement exists for the two programs with the 
important exception of at low speeds (20 knots and below) 

and M . At 40 knots airspeed, the two programs can be con- 
sidered to show reasonable correlation (although M differs 

W 

by a factor of six) and thus at this speed the C-81 model is • 
used throughout this study to investigate the importance of 
rotor degrees of freedom on derivative identification. It is 
questionable whether the C-81 model is representative at speeds 
below 40 knots and thus the Ref. 1 models are used at these 
conditions. However, the importance of rotor degrees of freedom 
must be assessed in the post-flight identification processing 
before pursuing identif icaion with the six degree of freedom 
Ref. 1 derivative model. 

2.3.2 CH-47 Production SAS Models 

The CH-47 SAS used throughout this investigation was supplied 
by LRC and is discussed in Ref. 12. The form of this SAS was sup- 
plied as a transfer function feedback which requires conversion to 
state variable representation. The procedure for accomplishing this 
and the resulting state variable representation is presented. 
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The SAS used throughout this investigation in transfer 
function form is presented in Eqs . (,2.13) through (2.18). 

Pitch SAS 


LONG = K 17.9s 


1. 98s+l 


Roll SAS 


Yaw SAS 


q U7.9s+1) * (. 364s+l) (4. 22s+l) » 

Kq = -.439 m/sec 

(17.3 in/sec) 


(2 


’LAT 


K 

047s + l) ’ K p = -* 127 m/sec 

(-5.0 in/sec) 


(2 


*RUD Tr 3.2s 1 

~F~ r (3725+1) * ( . 098s+l) ’ 


Velocity > 40 kts K. 


Roll Into Yaw 


(2 

.262 m/sec 
(-10.32 in/sec) 


-f 12 " K rp (4.Tr 5 H.i 7 > x rp =16 - 4 cn/xad/sec C2 

|(6 . 4 7 in/rad/sec) 

Sideslip SAS (converted to lateral velocity at 40 kts) 


-it 2 = tK v r r?W • A * • 00272 m/m/sec 

(.0326 in/fps) 


(2 


DCPT SAS (NOTE: Not implemented on LRC CH-47) 


LONG 


= +.0111 


(2 
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.13) 


.14) 


.15) 


.161 


. 17 ) 


.18) 


u 



Conversion of the transfer functions Eqs . (2.13) through (2.18) 
to state variable form is complicated by the presence of zeros 
in the numerator. Standard techniques are available and appear 
in many references (c.f. Ref. 16) for this type of conversion 
(i.e. with the presence of zeros in the numerator). This conver- 
sion is not unique and depends primarily on the type of represent- 
ation desired. The standard canonical form was selected for the 
structure of the transformation and is of the form shown in 
Eq. (2.19). 



y = [1 0 0 ] v (2.19) 

where 

n = order of denominator 

m = order of numerator 

v. r . . - . = canonical states 

y = output vector of transfer function state mode 

Complete details of this technique can be found in Ref. 16 (pp. 
182-190) . 

The state variable form for the CH-47 SAS presented in 
Eqs. (. 2 . 13 ) through (2.18) is shown in Eqs. (2.20) through (2.24). 
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Pitch SAS 


~v 


0 

1 

• 0 


V 1 


_ 1 . 2 8 9 

; 2 

= 

0 

0 

1 


V 2 

+ 

-3.2688 

v 3 


- . 0364 

- . 8179 

-3.041 


V 3 


8.8863 






— 



- 

_ 


Roll SAS 


5 L0NG V 1 

v, = -21.28 v. + 21.28 K *p 
4 4 p r 


6 LAT " v 4 

Yaw SAS 



0 

-3.189 


1 


v s" 

10.517 


_ V 6_ 


10.2 

-107.29 



8 nTTT , = V r 

RUD 5 
Roll Into Yaw 

v- » - . 2435 v 7 + . 2435 K • p 
7 7 rp v 

6 RUD = v 7 

Sideslip SAS (Lateral Velocity) 

v o = -5.75 v Q + 5.75 K *v 
8 8 v 

6 RUD = v 8 

DCPT SAS (V = 40 kts) 

S L0NG = + ' 0111 u 
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q 

( 2 . 20 ) 


( 2 . 21 ) 


( 2 , 22 ) 


(2.23) 


(2.24) 


(2.25) 



The complete closed-loop model with SAS can be written in 
state variable form as shown in Eq. (2.26). 


• 

X 


F GH C 


X 


G 


— 




+ 


V 


Vc A C _ 


V 


0 


where x represents the open-loop state vector of the appropriate 
model (10 DOF or 6 DOF) and v represents the state vector of the 
SAS (8 state variables) . 

Appendices A and B present the model structure of the com- 
bined state vector models augmented by the SAS in state variable 
form (i.e., Eq. (2.26)) used throughout this study. 

2.3.3 Input Control Design Models 

Two basic approaches available for input design at this time 
are based on linear models. These are: 

(1) Utilization of specialized algorithms which extremize 
the information on the parameters of interest (Ref. 8) ; 

(2) Eigenstructure (e.g. eigenvalues and eigenvectors) of 
the system; 

The model analysis of the previous sections presents a chal- 
lenging problem for input design because of the high order of the 
coupled dynamics models and the presence of significant effects 
of nonlinear aerodynamic and dynamic interactions. The approach 
to this CH-47 input design was therefore accomplished by the 
following method. 

(1) Evaluation of the inputs against the closed- 

loop eigensvstem corresponding to the 10 DOF (SAS-on) 
system to determine frequencies and phases of the 
inputs ; 

(2) Simulation evaluation of the inputs on 10 DOF simulated 
data to further optimize the identification of the 
desired parameters and to evaluate the magnitude of 
vehicle responses (which should not be "excessive”) . 
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This approach concentrates on input designs which will 
enhance the uncoupled 3 DOF model parameters, and is evaluated 
upon 3 DOF and 10 DOF (coupled) responses. It should be noted 
that this approach does not imply that cross - coupled derivatives 
are not identifiable from the subsequent flight test data, nor 
should only 3 DOF models be used for post-flight processing. 

Based upon the discussion above, the procedure employed for 
parameter identification is outlined in Figure 2.12. 



I 

1 ► 3 DOF DERIVATIVES 

(COMPARE WITH quasi-static) 


Figure 2.12 Simulation Procedure Employed for Identification 
of Uncoupled 3 DOF CH-47 Derivatives 

The 3 DOF model structure used in the identification pro- 
cedure including the appropriate SAS is presented in Eq. (2.27) for 
the longitudinal case and the baseline measurement equation is 
shown in Eq. (2.28). 
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L'.J 


The lateral model structure is shown in Eq. [2.29) and the 
baseline measurement equation in Eq. (2.30). 
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where 




(2.30) 


The longitudinal model includes the complete pitch SAS, 
whereas the lateral SAS, which contains 5 state variables, was 
reduced to the 2 state -variable SAS as shown in Eq. (2.29). This 
is justified since the frequency associated with the SAS state 
variables which are deleted is very high. The three SAS variables 
deleted are roll rate feedback (v^) » yaw rate feedback (v^.) » and 
sideslip ( v g)> with time constants of 0.047, 0.098, and 0.174 
respectively. It should be noted that the roll and yaw SAS time 
constants are in the vicinity of the rotor modes and, since the 
rotor is not modeled in this analysis, these SAS state variables 
are deleted. 

The effect of the reduced SAS in the lateral equations (Eq. 
(2.29)) is to modify the open-loop derivatives as indicated by the 
primed derivative values (e.g., Y" = Y +K Y6 DIin ) . The primed 
values are the parameters to be identified and the open- loop 
values are then derived (e.g., Y = Y - K Y5 min ) . 

V V V KUU 

2.4 INSTRUMENTATION AND DATA ACQUISITION 


Instrumentation and the aircraft data recording facilities 
constitute one of the factors influencing parameter identification 
accuracy* An extensive analysis of this factor is documented in 
Ref. 3, with further investigations published in Refs. 18 and 19. 
The modeling of the instrumentation of the VALT Research Aircraft 
and LRC data recording facilities follows the procedures and 
conclusions presented in Ref. 3. 
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Figure 2.13 presents an overall diagram of the instrumenta- 
tion modeling procedure used in this investigation. The CH-47 is 
shown with SAS operative, and the vehicle sensors and LRC data 
recording system models are shown in this figure. The models indi 
cated in Figure 2.13 include dynamic lags, scale factor matrices, 
random error, and bias errors on both system measurements and 
control measurements. 

The baseline values used to represent the sources of error 
indicated in Figure 2.13 are presented in Table 2.4. These values 
were obtained by examination of the individual error specifica- 
tion sheets associated with the instrumentation of the CH-47 VALT 
Program Research Aircraft, the PADS (Piloted Aircraft Data System) 
and the RAGS (Research Aircraft Ground Station Data Reduction 
System) data systems. Although these specification sheets present 
essentially overall system upper bound accuracy limits which are 
not representative of Gaussian standard deviations (lo) , the 
values presented in Table 2.4 were selected based upon best esti- 
mates which take into account those considerations addressed in 
Ref. 3. 

In addition to the baseline (nominal) error magnitudes 
presented in Table 2.4, both larger and smaller error magnitudes 
were investigated in this study. These error magnitudes will be 
further discussed and presented in the section which gives 
results of the simulation study CSection 3.1.1). 

2.5 FLIGHT ENVIRONMENT 


The flight test environment constitutes an important factor 
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tives. Its importance to helicopter identification is a result 
of low speed, low altitude flight and the high level of vehicle 
sensitivity to wind gust due to the speed dependent stability 
derivatives (e.g., M u and M^) . The flight test environment 
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Table 2,4 Nominal CH-47 VALT Flight Instrumentation Measurement Errors 
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important to identification accuracy can be considered to consist 
of three elements: random turbulence, wind shear, and ground 

effects. Each of these are discussed in the following. 

Random Turbulence . The random turbulence model used in this 
study is based upon theoretical and experimental information. The 
nature of random turbulence in the general case is such that 
present theoretical and experimental analysis techniques only 
approximate the actual random characteristics. The simplest model 
for wind gust is the well known Dryden model of atmospheric 
turbulence. Extensions and developments are continually being 
researched to improve upon this characterization [20] . Among 
those included are modeling effects directed toward the non- 
Gaussian nature of turbulence (in particular, "patchiness"), 
spatial and time variability, and nonisotropic representation. 
Although each of these modeling areas represents an improvement in 
the actual characterization of random turbulence, it is not the 
intent of this investigation to model wind gust, but to assess 
the influence of random wind disturbances on parameter identifi- 
j cation accuracy. For this reason, the Dryden model is chosen. 

; The following three assumptions are made in regard to the wind 
| turbulence model used in this investigation: 

• Isotropic (i.e., equal scale lengths and zero cross- 
spectra) 

• The turbulence field is "frozen" (i.e., time and spatial 
frequencies are related by the relative tip airspeed of 
the rotor) 

• The power spectral density for the horizontal velocity 
component is used for the lateral and vertical compo- 
nents 

The last assumption is believed to have no consequence on 
the general conclusions of this study. The vertical and lateral 
Dryden model spectrums differ slightly in comparison to the 
longitudinal axis, thus modifying the distribution of energy 
in the low frequency region. 



Equation (2.31) presents the Dryden gust spectral representa- 
tion (Ref. 20) . 


where , 


<i> u (w) 


u 


0) 7T 
0 


2 

[1+ (u/w 0 ) 2 ] 


* *v - <!*„ C2.31) 

assumed for this study 


u) = (fiR + V) / JL 

Q = rotor rotational speed, not the spatial frequency 


The turbulence scale length, l, is set to 122 meters. The 
spectral break frequency io is a function of airspeed, rotor tip 
speed, and correlation distance. (Rotor tip speed is used to 
modify airspeed since at hover, the break frequency goes to zero 
and, as computed in Ref. 12, gives an improved representation 
based upon simulation analysis.) 

Equation (2.32) presents the time domain version of the longi- 
tudinal gust model and is used throughout this investigation for 
modeling all three axes as discussed above. 


u = -u u + d) ?(t) (2.32) 

g O g o 

At hover, = 1.8 rad/sec and £(t) represents random white’ 
Gaussian noise with standard deviations selected such that the 
output of Eq. (2.32) has standardf deviations of o u# 

Table 2.5 summarizes definitions of the turbulence level 
assumed for this investigation. 


Mean Wind Profile (Wind Shear) . In addition to random tur- 
bulence, parameter identification accuracy may also be affected by 
the presence of a wind shear. This is particularly true in a 
landing approach where the wind shear effect is most pronounced. 

As the vehicle descends in the presence of a wind shear, the 
wind velocity decreases and this changes the relationship between 
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Table 2.5 

Definitions of Turbulence Level 



the inertial and airm'ass relative velocity vectors. Depending 
upon how a pilot handles the helicopter in the wind shear (i.e., 
does he hold 0, v^,, h?) and the vehicle's inertial properties, the 
measured accelerations will vary due to the wind shear. Parameter 
identification accuracy is, thus, dependent upon knowing the 
relationship between the inertial and relative velocity vectors. 

For this study, a detailed quantitative appraisal of wind 
shear effects on parameter identification accuracy is not made. 
Because of this, the proposed flight test plan is premised on 
testing only when significant shears are not present. 

Ground Effect . Ground effect is usually modeled as a 
correction to thrust coefficient as a function of hub altitude 
from the ground and hub velocity. This correction is typically 
insignificant above altitudes one-half to three- fourths of the 
rotor diameter. Above approximately 15 meter altitude, the CH-47 
can be expected to be out -of - ground effect. The procedure recom- 
mended for determining the influence of ground effect is to test 
the flight vehicle at two conditions -- in-ground effect (below 50 
meters) and out-of -ground effect (above 15 meters) . Parameters 
identified at each condition can be compared and differences in 
parameter values attributed to the proximity of the ground. 
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2.6 FLIGHT TEST INPUT DESIGN CONSIDERATIONS 


The design of flight test control inputs is complicated for 
rotorcraft due to the many constraints imposed upon the design. 
These constraints include simplicity for practical implementation, 
good excitation of modes for good identif iability , frequency con- 
straints to avoid excessive rotor excitation, state variable 
constraints to retain linearity, and pilot stabilization at 
various flight test conditions. In addition, the presence of 
a dynamic state variable SAS required for additional stability 
complicates any analytic design procedure since it is desired to 
identify open-loop derivatives with the SAS operative. Figure 2.14 
summarizes those considerations important to rotorcraft input 
design for parameter identification including the various con- 
straints imposed for the CH-47. The nominal flight test input 
design used in this effort is of the shape shown in this figure. 
Details of the actual inputs recommended will be presented in the 
next chapter. 
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Figure 2.14 Flight Test Input Design 




III. IDENTIFICATION RESULTS USING CH-47 SIMULATION 


This chapter presents the simulation results upon which the 
CH-47 flight test design is based. All simulation runs include 
the CH-47 production SAS with the roll rate feedback gain divided 
by a factor of two (except where noted) to eliminate fuselage/ 
rotor coupling as discussed in detail in Chapter II (Section 2.2.2). 
The procedure employed in the simulation studies follows the 
Technical Approach outlined in Figure 1.2. 

Table 3.1 provides a summary of the simulation conditions 
investigated in this study. The models employed at the 40 kt 
and 60 kt conditions are derived from the C-81 simulation pro- 
gram and the models used for the hover condition were supplied 
by LRC. The conditions indicated by the check (/) represent 
actual computer simulations. The previously mentioned computa- 
tional problems with C-81 precluded a detailed evaluation of 
all test points. 

3.1 LONGITUDINAL IDENTIFICATION RESULTS AT THE 40 KT CONDITION 

This section presents the longitudinal identification 
results at the 40 kt condition. This analysis is generated from 
the 10 DOF, 6 DOF, or 3 DOF models shown in Appendix A; and 
control inputs consist of either a longitudinal input (differ- 
ential collective) applied separately or in combination with a 
collective input. This section summarizes the various control 
inputs investigated and presents the final input designs recom- 

m A -Pr\ -r» -Pi i rrVi •+■ f a C + 
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Table 3.1 

Simulation Conditions Examined Leading to Flight Test Specification 
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Table 3.1 
(Concluded) 
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3.1.1 Control Input Design Procedure and Simulation Results 


The 3 DOF closed-loop model representing the CH-47 at the 
40 kt condition was used for designing inputs. The characteris- 
tic roots of the closed-loop system, natural frequency, damping, 
and period are presented below. 


Table 3.2 

CH-47 Closed-Loop Characteristic Roots (40 Kts) 



EIGENVALUE 

to 

n 

C 

T 

(SECS) 

Pitch Damping 
+ SAS: 

-- I. 977+j2.80 

3.43 

.572 

1.83 

Vertical/Long. : 

-. 517+j . 09348 

.525 

.984 

12.0 

Long. Divergent 
+ SAS: 

-. 000092+ j. 0859 

. 0859 

.0107 

73.0 

SAS: 

-.0468 

0 

— 

— 


The longitudinal input shown in Eq. C3.1) was initially se- 
lected such that the three frequencies were chosen to be those of 
the closed-loop natural frequencies. 

‘■'LONG = l®* 2 s ^ n 3.43t + 2.5 sin .525t + 2.5 sin .0859t (cm) (3.1) 


Application of this input to the closed-loop CH-47 model resulted 
in state variable responses of very large magnitudes; considerably 
greater than the small perturbations of the linear model. After 
severax simulation attempts with different centre! inputs, the 
input design shown in Eq. (3.2) was found to produce acceptable 
state response: 


^ LONG = A 1 sin 3 - 43t + A 2 sin ,424t C 01 ®) 


(3.2) 
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where, based on simulation in response evaluation, 


(2.5 

t 

< 

7.32 sec 

10. 0 

t 

> 

7.32 

(.64 

t 

< 

12.0 

I o 

t 

> 

12.0 sec 


The input used in Eq. (3.2) is such that four complete 
cycles occur at the frequency 3.43 rad/sec, and one complete 
cycle at the frequency 0.525 rad/sec. This input excites all of 
the longitudinal state variables with good amplitude, yet the per- 
turbations are within the linear model assumptions. 

The input design of Eq. (3.2) was used to generate data 
for use in the identification phase. The data were contaminated 
with random noise and the length of data used was 15 sec, sampled 
at 0.04 sec increments. 


The longitudinal model of Eqs . (2.27) and (2.28) was used 
to assess the effects of random noise on derivative prediction 
accuracy. This investigation also determines if it is possible 
to obtain accurate identification using only the measurements of 
Eq. (2.28). Random noise was added to the measurements and was 
chosen to be that used in the previous instrumentation require- 
ments study of Ref. 3. The random noise was selected to be 
the large noise case and is shown in Eq. (3.3). 

a = .0017 rad/sec 1 ' 
a „ = .0034 rad 

0 (XX 

v. ~ ~ 

a.. = .01 g's 

Z 


a 


A X 


.01 g's 
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The data base generated using the 3 DOF model shown in Eqs . 
(2.27) and (2.28) was used to identify the open-loop stability 
and control derivatives. The level of measurement noise used is 
as shown in Eq. (3.3) and the control input used is shown in 
Eq. (3.2). The maximum likelihood identification program (SCIDNT) 
was modified to accommodate the model structure shown in Eqs. 
(2.27) and (2.28).. 

Table 3.3 summarizes the results of the identification and 
presents the parameter start-up used, maximum likelihood esti- 
mated parameters and a la band, and the actual error. 


Table 3.3 

Results of Longitudinal 3 DOF Maximum 
Likelihood Identification 
(V = 40 Kts, High Noise Case) 


IDENTIFIED 

PARAMETER 

ASSUMED 

TRUE 

PARAMETER 

VALUE 

INITIAL 

GUESS 

ML 

ESTIMATED 

PARAMETER 

ML EST. 

1 <7 

(ABS) 

ACTUAL 

PAR. 

ERROR 

ACTUAL 

% 

ERROR 

X u 

• .0087 

‘-.01 

-.00399 

.00244 

.00742 

54.1 

Z u 

" 117 

.0 

-.1148 

.00805 

.0096 

1.9 

H u 

-.00705 

.0 

.00761 

.0001716 

.00056 

7.9 

*w 

.06069 

.0 

.09664 

. 01 308 

.03595 

59.0 

Z w 

-.6566 

-.5 

-.6248 

.01949 

.0318 

4.84 

H w 

.0121 

.003 

.01236 

.000814 

.00020 

1.62 

M q 

-1.329 

-1.0 

-1,3225 

. 006045 

.0065 

.49 

X ' 5 LQNG 

2.9484 

.0 

2.14896 

.45276 

.79944 

27.1 

Z,5 LQNG 

8.5380 

2.4 

8.4324 

.57684 

.1056 

1.24 

M,5 L0HG 

-15.5906 

-11.811 

-1.29941 

-15.5929134 

.02362 

.0151 




Examination of the actual error in the estimated deriva- 
tives of Table 3.3 reveals that all derivatives are identified 

with less than 101 error with the exception of X , X , and 

u w. 7 

X^long* These derivatives all occur in the X-force equation. The 
longitudinal acceleration A x provides most of the information 
regarding these derivatives and the level of random noise used 
on the A x measurement is nearly as large as the signal level. 
This accounts for the poor accuracy on the X-force derivatives. 
Also, it should be noted that the level of noise chosen for the 
A measurement (0.01 g's) represents the high noise case. Reduc- 
tion of this noise level to be more representative of light data 
can improve the accuracy further. 

In developing a control input for actual implementation, 
it is desirable to have simplicity in the shape of the input, 
yet retain good accuracy in identified parameters. The effect 
of simplifying the shape of the input shown in Eq. (3.2) is pre- 
sented below, along with the effect of noise level, lateral-to- 
longitudinal coupling, and rotor degrees of freedom. 

The time history response data is generated from the 3 DOF, 

6 DOF, and 10 DOF linear C-81 derivative model at the 40 kt 
condition. In each case, the complete CH-47 SAS is included. 

Table 3.4 summarizes the various sets of data generated for use 


Table 3.4 

Summary of Time History Data Generated for Use with SCIDNT 

(CH-47C C-81 Models, 40 Kts) 


MODEL USED 
IN DATA 

rPMC n ATTAM 
OuiLlVi 1 lull 

(SAS-0N) 

MULTIPLE CYCLE INPUT 

1 CYCLE INPUT 

mcai tptoi r 

IILOLlCiiUUL 

NOISE 

M amt ma i 
liUru Unu 

NOISE 

U T All 
HiOil 

NOISE 

LOW NOISE 

3 DOF 

/ 

/ 

/ 

/ 

6 DOF 



/ 


10 DOF 


/ 

/ 

/ 
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in the parameter identification program (SCIDNT) . Table 3.5 
summarizes the noise levels used for the cases described in 
Table 3.4. 


Table 3 . 5 

Noise Levels Used in Data Simulation 


TIME HISTORY 
MEASUREMENT 

NOISE LEVEL 

STANDARD DEVIATIONS 

NEGLIGIBLE 1 

NOMINAL 2 

HIGH 

q rad/sec 

.0000173 

.00054 

.00173 

0 rad 

.0000173 

.0011 

.0034 

A x g's 

.0001 

.0033 

,.01 

A 2 9 ' S 

.0001 

.0033 

! .01 


Negligible noise = o High * 100 
Nominal noise = o^.. ^ * 3 


Two different longitudinal control inputs (differential 
collective) were investigated: (1) a multiple cycle input, and 

(2) a single cycle input. These inputs are shown below. 

Multiple Cycle Input 


4 cycles 


1 cycle 


'LONG 


= A, sin 3.43t + A~ sin .525t 


cm 


(3.4) 


(2.54 

t 

< 

7.32 

sec 

(0.0 

t 

> 

7.32 

sec 

(0.635 

t 

< 

12.0 

sec 

( 0.0 

t 

> 

12.0 

sec 
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1 cycle 


Single Cycle Input 


5 L0NG 


1 cycle 


sin 3. 43t 


+ 



A 2 sin . 525t 


cm 


(3.5) 


_ 12.54 

. t 

< 

1.8 sec 

( 0.0 

t 

> 

1.8 sec 

(0.635 

t 

< 

12.0 sec 

(070 " 

t 

:> 

12.0 sec 


These inputs contain two frequencies (to.. = 3.43, = 

.525 rad/sec) which were chosen to be equal to the natural fre- 
quencies of the closed-loop system. The multiple frequency 
input of Eq . (3.4) contains 4 complete cycles of high frequency, 
which is designed to excite the short period mode of the closed- 
loop system. The low frequency input contains 1 cycle and is 
designed to excite the longitudinal/vertical coupled mode and 
long period mode. 

The single cycle input shown in Eq. (3.5) is used in this 
study to assess whether accuracy in derivative identification is 
impaired as compared to the multiple cycle input. The single 
cycle input is preferred in a flight test program due to its 
simplicity (essentially a doublet) and if accuracy is not 
degraded, is the recommended input for implementation. 

SCIDNT was used to assess the effects on derivative identifi- 
cation accuracy of: (1) noise level, (2) longitudinal-to-lateral 

coupling, (3) rotor degrees-of-freedom, and (4) control input 
simplification. In each case, the model structure programmed in 
SCIDNT contains only the longitudinal 3 DOF derivative model with 
longitudinal SAS. The form of the equation used in SCIDNT is 
presented in Eqs . (2.27) and (2.28). The influence of: (1) noise 
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level, (2) longitudinal- to-lateral coupling , (3) rotor degrees-of- 
freedom, and control input simplification, are determined. (Refer 
to Table 3.4 for a summary of the cases investigated.) 

Effect of Noise Level on Derivative Identification Accuracy 

The effect of noise level on derivative identification ac- 
curacy is summarized in Table 3.6. Data are generated from the 
3 DOF longitudinal model for the three cases — negligible noise, 
nominal noise, and high noise. For each case, the identified 
derivative values and error in the estimate are presented. As 
shown in Table 3.6, all derivatives are identified with less than 
10% error (excellent accuracy) for the nominal noise case. The 
high noise case causes further degradation in the accuracy of the 
derivatives X , X w and The reason for this is that the 

term (g cos 6 q )0 contributes significantly to the u time history; 
thus, the derivatives X , X w , and ^6 l 0NG are difficult t0 iden- 
tify when noise is high. The results of Table 3.6 indicate that 
for the noise levels typical of flight data (nominal noise case), 
accuracy will be excellent. 

Effect of Longitudinal- to-Lateral Coupling on Derivative 

Identification Accuracy 

The effects of lateral coupling on longitudinal identifica- 
tion accuracy were determined by generating time history data 
using the 6 DOF coupled longitudinal and lateral model (with 
complete SAS) and attempting to identify only the 3 DOF longi- 
tudinal derivatives. Again, the 3 DOF (with SAS) model was used 
in SCIDNT . Table 3.7 summarizes the identification results for 
the high noise case and shows identified derivatives from data 
generated by the 3 DOF simulation and the 6 DOF simulation. 

A comparison of 3 DOF and 6 DOF results (Table 3.7) reveals 
that lateral coupling altered the prediction accuracy in X u 
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Table 3.6 

hffect of Noise Level on Derivative Identification Accuracy 
(Data Generation' from 3 DOF Model with SAS, CH-47C, 40 KTS) 
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NOMINAL NOISE 




Table 3.7 

Effect of Lateral Coupling and Rotor DOF on Longitudinal 
Derivative Identification Accuracy (High Noise Case) 
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Errors Due to Errors Due to Errors Due to 

Random Noise Random Noise + Random Noise + 

Lat. Coupling Lat. Coupling + 

Rotor OOF 




and M w and that the accuracy in other derivatives is essen- 
tially the same. The derivative M w is affected by coupling 
because its contribution to the q time history is very small 
at this condition (40 kts) using the C-81 derivative model. 

(Note that the C-81 value for M w is .0037, whereas the Ref. 1 
derivative for M w is .0207; five times as large.) The finding 
that lateral coupling has a negligible effect on longitudinal 
derivative identification is due to two factors: (1) longi- 

tudinal inputs do not excite the lateral - directional degrees of 
freedom, and (2) lateral -direct ional motions are suppressed by 
the stability augmentation system. 

Effect of Rotor Degrees- of- Freedom on Derivative Identi - 
fication Accuracy 

The effect of rotor DOF on derivative identification ac- 
curacy was determined by generating time history data using the 
10 DOF fuselage-flap derivative model (with complete SAS) and 
identifying only the 3 DOF longitudinal derivatives. Table 3.7 
summarizes the identification results for the high noise case 
and shows the results when using 3 DOF, 6 DOF and 10 DOF simula- 
tion data (SAS-on). Notice that there is further degradation 
of ML, (61.1% error), M„ is now at 10.1% error, and the X- 
force derivative (X w and x Slong) errors are a combined result 
of high noise level, lateral coupling, and rotor DOF. 

Table 3.8 summarizes the effect of both lateral coupling 
and rotor degrees-of- freedom on derivative identification accur- 
acy for the nominal noise case . Identification results from 
the 10 DOF simulation shows the degradation due to lateral coupl- 
ing and rotor DOF. Again, M w is identified to 48.7% error; 
part of this error due to lateral coupling and part due to 
rotor DOF. The X-force derivatives X w and ^SlqnG are ac * en “ 
tified to 36.7% and 20.8% error, respectively, and these errors 
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Table 3.8 


Effect of Lateral Coupline and Rotor DOF on Derivative 
Accuracy (Nominal Noise Case) 


PARAMETER 

C-81 

TRUE 

PARAMETER 

VALUE 

DATA FROM 3 DOF 
SIMULATION 

0ATA FROM 10 OOF 
SIMULATION 

MAXIMUM 

LIKELIHOOD 

ESTIMATED 

PARAMETER 

% ERROR 

MAXIMUM 

LIKELIHOOD 

ESTIMATED 

PARAMETER 

% ERROR 

X u 

-.01141 

-.01078 

5.3 

-.0121 

6.7 

Z u 

-.1244 

-.1182 

5.0 

-.1203 

3.3 

H u 

.00705 

.00735 

4.7 

.00676 

4.2 

X w 

.06069 

.0531 

12.5 

.03837 

36.7 

z w 

-.6566 1 

-.6366 

3.0 

-.6356 

3.2 

\ 

.0121 

.0112 

8.1 

.00630 

48.7 

% 

-1.329 

-1.318 

1.0 

-1.202 

9.6 

Xfi Long 

2.9484 

3.2508 

10.3 

3.5592 

20.8 

Zfi Long 

8.5380 

8.2656 

3.2 

7.9392 

7.0 

^Long 

-15.5906 

rl5.5748 

0.1 

-15.2953 

1.9 


Errors Due to Errors Due to 

Random Noise Randan Noise + 

Lat. Coupling + 
Rotor OOF 


result because of the small influence these derivatives have in 


V ~ V“ ~ ~ 1 

Ai.O J\ xO i CC w\|UCk w aWu < 


T n t?i imm n T*Ar 

X A* O Miltil 


r\ro^i rtinn +■ *r c\ q for thp 

j y ±- — s ~ *- ’ — w — “ w ** — 


3 DOF and 10 DOF results are similar for both noise levels inves- 


tigated. 
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Effect of Control Input Simplification on Derivative Iden - 
tification Accuracy 

Two different control inputs were investigated for the dif- 
ferential collective control C^lonG^ : ^ a mu ltiple cycle 

input, and (2) a single cycle input. The precise form of the 
inputs used are shown in Eqs. (3.4) and (3.5). Table 3.9 summar- 
izes the identification results and accuracies obtained for the 
two inputs using both the 3 DOF simulation and the 10 DOF simu- 
lation. The results using the 3 DOF simulation data indicate 
that the single cycle input yields results nearly as good as 
the multiple cycle (4 cycles) input. The results using the 10 
DOF simulation also indicate that a single cycle input is satis- 
factory. 

Table 3.10 summarizes the eigenvalues from the identified 
derivatives for the nominal noise case with the single cycle 
input. Eigenvalues are shown from derivatives identified using 
the 3 DOF simulation and the 10 DOF simulation and compared with 
the Ref. 1 3 DOF values. Table 3.10 also shows that the eigen- 
values from the identified derivatives using the 3 DOF and 10 DOF 
simulations are similar and that their agreement with the "true" 
values is excellent. 


An important element of the overall identification is the 
standard deviations for parameters and measurement noise. Esti- 
mated parameters from flight data are themselves incomplete since 
their usefulness is limited without accuracy confidence bands. 
Table 3.11 shows the estimated la confidence bands, as determined 
via the maximum likelihood information matrix. In all cases, 


except M of 
either la or 


xable 3.11, the actual errors oic i ©presented by 
2a bands. A 2 a band is an upper bound on the errors. 


Table 3.12 shows 3a confidence bands for the nominal case 
where derivatives were identified using the 10 DOF simulation data. 
The effects of lateral coupling and rotor DOF degrade the accuracy 
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Table 3.9 

Effect of Control Input Simplification on Derivative Identification Accuracy 
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Table 3.10 


Eigenvalues from Identified Derivatives, Nominal Noise, 

1 Cycle Input 



C-81 

TRUE EIGENVALUES 

USING DERIVATIVES 
IDENTIFIED FROM 
3 DOF SIMULATION 

USING DERIVATIVES 
IDENTIFIED FROM 
10 DOF SIMULATION 

PITCH DAMPING 

-1.977 + j 2.804 

-1.974 + j 2.806 

-1.904 + j 2.775 

VE RT I CAL/ LONG ITU0I NAL 

-.5174 + j .0935 

-.5101 + j .0966 

-.524 + j .0451 

LONGITUDINAL DIVERGENCE 

-.00092 +j .0859 

- . 00027 + j . 0904 

.00063 + j .0929 

SAS 

-.0468 

-.0468 

-.0473 


Table 3.11 

Parameter Standard Deviations for the Nominal Noise Case, 
1 Cycle Input (Derivatives Identified from 
; 3 DOF Simulation _Data, 40 Knots) 


PARAMETER 

C-81 

TRUE 

PARAMETER 

VALUE 

MAX. LIKELIHOOD ESTIMATED 

ACTUAL 

ERROR 

Tabs. | 

PARAMETER 

ESTIMATE 

la* 

X u 

-.01141 

-.01009 

.0009 

.0013 

Z u 

-.1244 

-.1155 

.0034 

.0089 

"u 

.00705 

.00745 

.000082 

.00039 

*w 

.06069 

.05970 

.0057 

.001 

Z w 

-.6566 

-.5501 

.0083 

.0165 

"w 

.0121 

.01207 

.00043 

.00007 

M q 

-1.329 

-1.323 

.0053 

.0054 

^Long 

2.9484 

3.1524 

.276 

.204 

^Long 

8.5380 

8.2236 

.3204 

.312 

^Long 

-15.5906 

-15.6102 

.02480 

.01966 
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of the confidence bands. A 3a to 6a confidence band more accu- 
rately reflects the actual errors in the derivatives when using 
the 10 DOF simulated data. 


Table 3.12 

Parameter Standard Deviations for the Nominal Noise Case, 
1 Cycle Input (Derivatives Identified from 
the 10 DOF Simulation Data, 40 Knots) 



C-81 

TRUE 

MAX. LIKELIHOOD ESTIMATED 

ACTUAL 

PARAMETER 

PARAMETER 

VALUE 

PARAMETER 

ESTIMATE 

3a* 

ERROR 
I ABS. | 

X u 

-.01141 

-.01809 

.003 

.0067 

Z u 

-.1244 

-.12108 

.013 

.0033 

M u 

.00705 

.00686 

.00020 

.00020 

*w 

.06069 

.05226 

.018 

.0084 


-.6566 

-.6378 ' 

.024 

.0187 

M w 

.0121 

.00650 

.0013 

.0056 

M q 

-1.329 

-1.2148 

.016 

.114 

^Long 

2.9484 

3.2796 

.864 

.3324 

^Long 

8.5380 

7.9368 

.936 

.6012 

^Long 

-15.5906 

-15.2913 

.0748 

.2992 


3a bands shown accurately represent only 6 out of 10 actual 
parameter errors 


The measurement noise covariance is estimated in the 
maximum likelihood program (SCIDNT) and the values for the nom- 
inal noise case, 1 cycle input are shown in Table 3.13. As 
shown in Table 3.13, those obtained using the 3 DOF simulation 
data and those obtained using the 10 DOF simulation are accur- 
ately identified. 
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Table 3.13 


Measurement Noise Standard Deviations Estimated From 
SCIDNT, Nominal Noise Case, 1 Cycle Input 


MEASUREMENT 

NOISE 

TRUE* 

STANDARD 

DEVIATION 

ESTIMATED STANDARD DEVIATION 

3 DOF 

SIMULATION 

CASE 

10 DOF 
SIMULATION 
CASE 

(7^ rad/sec 

.00058 

. 00073 

. 00082 

a Q rad 

.00115 

.00133 

.00134 

aA z g's 

.0033 

.0032 

.0032 

oA x 9' s 

,.0033 

i .0033 

; .0035 


* 

Used in the simulation 


3.1.2 Application of Model Structure Determination 

The procedure used throughout this contract is to treat the 
decoupled problem (i.e., longitudinal separate from lateral). The 
results of the studies to-date indicate that the presence of lateral 
coupling provides very little effect on longitudinal identification 
accuracy. However, the CH-47 flight test vehicle could possibly 
contain greater coupling than the analytic simulation models and, 
thus, it would be desirable to determine the degree of coupling in 
the flight data itself. Also, the question of identification of 
the cross-coupling derivatives is raised. 

Since one of the objectives of this effort is to provide 
stability and control derivatives for flight control systems 
design, one consideration is to accurately identify the longitudi- 
nal and lateral decoupled derivative arrays. In addition to this, 
it is required to: (1) determine from the flight data the degree 

of lateral coupling, and (2) identify the required cross -coupled 
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derivatives if coupling is found to be significant. The assessment 
of the amount of lateral coupling can be accomplished through the 
use of a model structure determination program. 


A model structure determination program was applied to the 
6 DOF model (with SAS) for a multiple cycle input case [i.e., Eq. 
(3.4)]. The model -structure determination technique operates on 
the state variables and accelerations of the system and produces 
the "best model" (i.e., simplest model) required to explain 
the data. This model building procedure is based upon a step- 
wise regression procedure which tests the significance of each 
possible derivative term, adds or deletes the terms from the 
model, and provides regression estimates of those derivatives 
found to be significant. The output of this procedure is a 
reduced order model which explains almost all of the variation 
occurring in the time history data (e.g., 991 of the variation). 


Application of the model structure determination program 
to time history data from the 6 DOF model (with the 8- state 
variable SAS) resulted in the model shown in Eq. (3.6). 
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(3.6) 


The derivative terms provided the following explanation of 

the time history data: 

REQUIRED % CONTRIBUTION 

DERIVATIVE TERMS to X, Z, or q 


X-Force: 

X , 
u’ 

X w> 

X q’ ^LONG 

98.82 

Z-Force : 

z , 

u’ 

Z W 

Z q’ Z6 L0NG 

99.95 

q-Equation:* 

M , 
u’ 

M q’ 

m§ long 

99.75 


*M was found to provide a negligible contribution to 
w the q time history. 
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The lateral terms are thus shown to provide a very small con- 
tribution to the longitudinal motion for the control input used. 
Notice that for the q-equation, the derivatives M u , M , and 
M<$long contribute 99.75% and contributes a very small 

amount. This is probably the reason that the identification 
accuracy of. M w deteriorates when lateral degrees-of- freedom 
and rotor degrees-of - freedom are present (see Tables 3.7 and 
3.8). Also, the M w in Eq . (3.6) is from C-81 and is one-fifth 
the size of the LRC derivatives discussed earlier. 

This example demonstrates a viable procedure for examina- 
tion of actual flight data to assess the degree of coupling in 
any time history response and the derivatives required for 
identification. 

Identification of Cross-Coupling Derivatives 

Although the simulation results clearly indicate that the 
SAS decouples the longitudinal and lateral response, it may be 
desirable to attempt identification of the cross-coupled deriva- 
tives. The procedure recommended for accomplishing this includes 
the following steps: 

(1) Generate one maneuver of data using longitudinal 
inputs . 

(2) Generate one maneuver of data. using lateral inputs. 

(3) Apply model structure determination program to deter- 
mine the significant derivatives for each maneuver. 

The longitudinal input (step 1) will provide information 
primarily for the longitudinal uncoupled derivatives. The 
lateral input (step 2) will provide information primarily for 
the lateral uncoupled derivatives. These derivative arrays 
can be identified using SCIDNT. Once these arrays are identi- 
fied, then the longitudinal input maneuvers (step !)• can be 



used with the lateral (uncoupled) derivatives identified from step 2 
and the longitudinal- to- lateral coupled derivatives identified. 
Likewise, the lateral -to -longitudinal (coupled) derivatives can be 
identified from the maneuver of step 2 using the longitudinal 
derivatives identified in step 1. 

This procedure reduces the problem of overparameterization, 
reduces the number of parameters required for each SCIDNT run, 
and should yield all the information required for flight con- 
trol system design. 

3.1.3 Final Control Input Design and Simulation Results 

The control input designs of the previous sections are 
based upon a differential collective input ( 5 lq N q) only. This 
section presents the results_of identification where control 
inputs for both longitudinal and collective stick are applied 
simultaneously in the maneuver. The results show improvement 
in derivative identification accuracy over the one input designs, 
and represent the final recommended inputs for the 40 kt longi- 
tudinal case. This section also presents the time history re- 
sponse for the 40 kt longitudinal case. 

Figure 3.1 shows the time history response due to a com- 
bined longitudinal C^loNG^ anc * collective (<$ C q LL ) control input. 

Each input contains two frequencies representing the short and 
long period modes of the closed-loop system as discussed previ- 
ously. The solid line response shown in Figure 3.1 is from the 
10 DOF rotor/fuselage/SAS model shown in Appendix A and the plus 
(+) symbols are from the 3 DOF decoupled longitudinal model 
(with' longitudinal SAS) also shown in Appendix A. The 3 DOF 
response is shown to accurately duplicate the 10 DOF response, 
thus indicating that a 3 DOF model would be sufficient for 
derivative identification purposes (rotor DOF and lateral coupl- 
ing are shown to have a negligible effect for the control inputs 
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10 OOF LINEAR MODEL 
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Figure 3.1 (Continued) 








CALCULATED 



Figure 3.1 (Concluded) 







employed). The lateral state variables have a negligible response 
demonstrating the effectiveness of the SAS in providing decoupled 
motion. 

Both the 10 DOF and 3 DOF simulation data of Figure 3.1 
were contaminated with the nominal random noise level and deriva- 
tives identified. The dual input identification results are shown 
in Table 3.14. The identified parameters using the 10 DOF simula- 
tion data show accuracy levels only slightly different compared 
with the results from the 3 DOF simulation data. These results are 
clearly superior to those obtained from the single input maneuvers 
^LONG^ were shown in Table 3.8. In particular, the 

derivative M w shows only 181 error compared to 48.7% error 
for the single input maneuver. 

Figure 3.2 shows the time history comparison of the 3 DOF 
identified derivative model of Table 3.14 and the 10 DOF simula- 
tion data with random noise contamination. Although most deriva- 
tives in Table 3.14 are identified to less than 10% error, M 

’ u 

and M w are identified with 18% error and are primarily respon- 
sible for the diverging fit shown occurring after 10 sec in 
Figure 3.2. This discrepancy in the time history match occurs 
because the long period root is nearly neutrally stable and 

M, and M have the most direct influence on this character- 
u w 

istic mode. As a result, the time response of the identified 
derivative model is very sensitive to small errors in M u and 
M.. . The eigenvalues for the simulation model, and identified 
derivative model using 3 DOF and 10 DOF simulation data, are 
shown in Table 3.15. 

3.2 LATERAL DIRECTIONAL IDENTIFICATION RESULTS AT 40 KTS 

Based upon the results obtained for the longitudinal 40 kt 
condition, the lateral-directional input design was chosen in a 
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Table 3.14 

Effect of Lateral Coupling and Rotor DOF on Derivative Accuracy — Two Control Inputs 

(40 KTS, Nominal Noise Case) 
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Table 3.15 


Eigenvalues of Identified Derivatives, 40 Knots 
Nominal Noise, Two Control Inputs 



C-81 

TRUE EIGENVALUES 

USING DERIVATIVES 
IDENTIFIED FROM 
3 DOF SIMULATION 

USING DERIVATIVES 
IDENTIFIED FROM 
10 DOF SIMULATION 

Pitch Damping 

-1.977 + j 2.804 

-1.969 + j 2.807 

-1.926 + j 2.779 

Vertlcal/Longl tudinal 

- -.'5174 + j .0935 

-.5166+j .0822 

-.526 + j .0691 

Longitudinal Divergence 

-.00092 + j .0859 

- . 00073_+j .0976 

-.00142+ j .1016 

SAS 

-.0468 

-.0478 

-.0483 


similar manner. A lateral C^lat^ an£ * a P ec * a l ^ ^ RUD ) i n P ut were 
chosen with two frequencies which occur at the short period and 
long period (Dutch roll] modes of motion. The amplitudes of the 
control inputs were chosen in an analogous manner to the longi- 
tudinal and were chosen such that the state variable response con- 
tains good excitation yet remains in a linear regime. Equations 
(3.7) and (3.8) present the control input selected for the 40 kt 
lateral-directional identification condition. 


S LAT = A 1 sin 3,5t + A 2 sin ' 319t ! (cm) ‘ (3.7) 

a- _ f3. 81 t < 1. 8 

1 ~ I 0 . 0 t > 1. 8 

A 2 = I. 25_ 


5 RUD 


B. 


sin 3.5t + 


B 2 sin 0.319t 


(cm) tes) 


B 


1 


B. 


{ 

{ 


-3.81 

0.0 

-.25 

0.0 


t > 6 and t < 7.6 

t < 6 and t > 7.8 

t > 6 

t < 6 


(3.8) 
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An earlier chapter discussed the problem of rotor- fuselage- 
SAS coupling as determined from the C-81 linear 10 DOF derivative 
model with the production SAS. This coupling was illustrated in 
Figure 2.8 which showed that both the rotor root locations and 
lateral short period root locations change when the CH-47 produc- 
tion SAS is included. It was also shown that this coupling 
could be essentially eliminated by reducing the CH-47 production 
SAS roll rate feedback gain by a factor of two. Table 3.16 
presents the eigenvalues of the various linear models derived 
from the C-81 nonlinear simulation program at the 40 knot condi- 
tion (SAS-off) . The effect of the CH-47 production SAS on the 
eigenvalues of the various models is shown in Table 3.17. 

Included in Table 3.17 are the eigenvalues with the roll rate 
feedback gain (K^) divided by a factor of two. This SAS modi- 
fication is shown to reduce the rotor coupling and place the 
lateral short period roots on the real axis as was illustrated 
in Figure 2.8. 

The lateral identification for the 40 knot condition in each 
case uses the 3 DOF lateral model with complex lateral SAS (5-state 
SAS), shown in Appendix A. The jnodel structure used in the identi- 
fication is as shown in Eq. (2.29) and contains two state variables 
for the SAS and assumes constant feedback gains for the remaining 
three states, as discussed in Chapter II. Identification results 
are presented for two cases: (1) simulation data employing the 

CH-47 production SAS and (2) simulation data generated with the 
CH-47 production SAS with the roll rate feedback gain K p divided 
by a factor of two. 

Table 3.18 presents the identified parameters for the 
lateral case where the data was generated from the 3 DOF lateral 
model with the 5-state variable CH-47 production SAS (Kp = -5.0). 
Those parameter values of Table 3.18 denoted with the prime 
represent the combined open- loop derivatives plus the constant 
feedback SAS contribution as shown in Eq. (2.29). The level of 
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Table 3.16 


Eigenvalues of the Various Linear Open-Loop 
Models Derived from C-81 (40 Kts) 


MODE 

3 DOF 
UNCOUPLED 

6 DOF 

10 DOF 

Pitch 

Damping 

-1.605 

-1.535 

-1.517 

Vert. /Long. 

-. 320+j .211 

-.310+j .261 

-.310 1 + J .261 

Long. 

Divergence 

.2512 

.2370 

.2377 

Roll 

Damping 

— 

-1.145 

-1.186 

Dutch Roll 

— 

+. 0905+ j .4198 

.0877+4.4159 

Spiral 

— 

-.04816 

- . 04828 

Flap 

-- - 

— 

- 18 . 08+ j 7.115 

Regressing 

— 

— 

-17.86+j7.256 

Flap 

— 

... 

-18. 134+ j 43. 53 

Advancing 

““ “ 

— 

-18.093+j43.44 


Table 3.17 

Eigenvalues of the Various Linear Models Derived from C-81 
with CH-47 Production SAS (40 Kts, Including DCPT Feedback] 


MODE 

3 DOF UNCOUPLED 
(PRODUCTION 
SAS) 

6 DOF 

(PRODUCTION 

SAS) 

10 DOF 
(PRODUCTION 
SAS) 

10 DOF 

K ■+ 2 
P 

Pitch Damping 

-1.977+42.804 

-1.974+42.817 

-1.914+j2.845 

-1.9138+42.837 

Vert. /tong. 

-.5174+j.09348 

-.5202+4.0923 

-.5199+4.0914 

-.527+4-0873 

Long. Divergence 

-.009214+4.0859 

-.000111+4.08517 

-.000183+4.0874 

.000866+4.0873 

Roll Damping 

-3.490+3.2385 

-3.491+4.2391 

-4.187+41.369 

-2.774, -1.156 

Dutch Roll 

- . 1158+_j . 3012 

-.11067+j.2998 

-.1122+4-2992 

-.2095+4.393 

Spiral 

-.00124 

-.00255 

-.00255 

-.0241 

Flap 

— 

— 

-9.66+42.34 

-15.7+45.43 

Regress'i ng 

— 

— 

- 18 . 7 + 45.49 

0 10 7C 



Flap 

— 



-18. 18+4 43. 53 

-18. 2+443. 57 

Advancing 



-18. 33+443. 61 

-18. 1+443.5 
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Table 3.18 


Identified Derivatives from the Lateral Simulation 
with CH-47 Production SAS (40 Kts , Nominal Noise) 


PARAMETER 

C-81 

TRUE 

PARAMETER 

VALUE 

MAXIMUM 

LIKELIHOOD 

ESTIMATED 

PARAMETER 

% ERROR 

V 

-.044544 

-.046323 

4.0 

V 

-.030272 

-.034268 

13.2 

n;* 

.013891 

.0139501 

.43 

v 

-2.074 

-2.3472 

13.2 

v 

-3.6416 

-4.0625 

11.6 

% 

-.37475 

-.57287 

52.9 


-.30727 

-.21590 

29.7 

L r 

.529015 

.73681 

39.3 


-1.99723 

-2.2939 

14.9 

Y<s LAT 

21.660 

14.9004 

31.2 

L(S lat* 

22.0827 

25.02244 

13.3 

n $lat 

2.10433 

3.151496 

49.8 

y<s rud 

-.66036 

-1.094640 

65.8 

l<s rud 

-2.3150 

-3.262008 

40.9 

N( Srud* 

7.40472 

8.82008 

19.1 










random noise contamination on the measurements of Eq. (2.30) 
represents the nominal case shown in Table 2.4. 

Figure 3.3 shows the time history responses of the identified 
derivative model of Table 3.18 and the simulated 3 DOF model. The 
time history match is shown to be an excellent fit even though a 
number of the identified derivatives (Table 3.18) are identified 
with large errors. Those derivatives of Table 3.18 denoted with 
an asterisk are considered primary and contribute most to the time 
response, thus resulting in good accuracy (<20$ error) . The remain- 
ing derivatives contribute little to the response and, thus, are 
difficult to identify. 

The lateral simulation model used. for generating the time 
history data was then modified by dividing the roll rate feedback 
gain (K ) by a factor of two. For this case, the time histories 
used in the identification were not contaminated by random noise to 
assess whether the derivative prediction accuracy could be improved. 
Table 3.19 shows large percent errors are still found for the 
secondary derivatives. Since zero noise contamination is used, 
these errors most likely are a result of using the approximate 
representation to the CH-47 SAS (i.e., the 2-state variable of 
Eq. (2.29) model rather than the 5-state variable model used in the 
simulation shown in Appendix A) . 

Table 3.20 presents the identified open-loop aerodynamic 
derivatives which are computed from Table 3.19 using Eq. (3.9). 
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Table 3.19 


Identified Derivatives from the Lateral Simulation 
with CH-47 Production SAS with Roll Rate Feedback 
Gain Divided by a Factor of Two (40 Kts) 

(No Measurement Noise) 


PARAMETER 

C-81 

TRUE 

PARAMETER 

VALUE 

MAXIMUM 

LIKELIHOOD 

ESTIMATED 

PARAMETER 

% ERROR 

N 

-.044544 

-.045790 

2.8 

L C 

-.030272 

-.0328068 

8.4 

»; 

.013891 

.0139692 

0.5 

% 

-1.2474 

-1.3207 

5.87 

s 

-2.2394 

-2.3370 

4.35 

% 

-.24113 

-.35453 

47.0 


-.30727 

-.24610 

19.9 

l ; 

.529015 

.68495 

29.46 

n ; 

-1.99723 

-2.2303 

11.7 

y<s lat 

21.660 

13.8720 

36.0 

l<s lat 

22.0827 

23.35472 

5.76 

n<s lat 

2.10433 

2.83925 

34.9 

y<s rud 

-.66036 

-.948276 

43.6 

L <SruD 

2.3150 

-2.993032 

29.3 

n< srud 

1 : 

7.40472 

i 

8.74173 

» - — 

18.0 







N = 
P 

% 

- K 

P 

N 5LAT 

(f) 

Y = 
r 

y ; 

- K 

r 

y srud 

(g) 

L r = 

L r 

- K r 

L 5RUD 

(h) 

N r " 

N r 

- K r 

N8 RUD 

(i) 

K v = 

.00272 : 

m/m/sec 



Kp = -.0635 m/sec 
K r = - . 262 m/sec 


Table 3.20 

Lateral Aerodynamic Derivatives Obtained 
from Equation (3.9) using Identified Values 
of Table -3.T19 (40 Knots) 
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Figure 3.4 Time History Comparison of the Identified Lateral Derivative Model 
with the Simulation Data (40 Kts , CH-47 Production SAS Gain K t 2 
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near perfect which further demonstrates that the derivatives of 
Table 3.19 identified with, large errors have a negligible influence 
on the response. The effect of reducing the CH-47 production SAS 
roll rate feedback gain and, thus, the rotor- to-airframe coupling 
is best seen by a comparison of the lateral acceleration responses 
(a ) in Figures 3.3 and 3.4. Reducing the gain improves the match 
between the 3 DOF simulation and the identified derivation model a 
responses . 

3.3 LONGITUDINAL IDENTIFICATION RESULTS AT THE HOVER CONDITION 

This section presents the results of the simulation studies 
for the longitudinal hover condition. The control input designs 
for the hover case employ single-axis inputs to constrain the 
appropriate state variables (either u or w) to retain linearity. 
This requirement was discussed earlier and is required due to 
the nonlinear variation at hover of the velocity derivatives 
M ' and M w , as shown i-n Figure 2.5. In addition to the above, 
the effect of data sample length and random wind gust on identi- 
fication accuracy is presented. 

3.3.1 Longitudinal Control Input Identification Results -- Zero 
Turbulence 


The identification procedure followed throughout the longi- 
tudinal hover simulations is as follows: (1) identify deriva- 

tives associated with longitudinal velocity, pitch rate, and 
longitudinal control, with excitation provided by a longitudinal 
control input ^LONG an ^ ^ identify the remaining vertical 
velocity and collective control derivatives, with excitation 
provided by a collective control input. Independent excitation of 
longitudinal velocity with vertical velocity fixed is required to 
avoid nonlinearity in the derivatives, as discussed in Section 2.2. 
This independent state excitation is possible at hover primarily 
because the SAS stabilizes the vehicle and essentially decouples 
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airspeed and pitch from vertical motion. Also, many of the cou- 
pling effects in the longitudinal equations are small at low speeds 
(particularly at hover and 2Q knots) . The degree to which this 
decoupling is achieved is best seen by examination of a time his- 
tory response. 


Figure 3.5 shows the time history response from the 6 DOF 
hover derivative model supplied by LRC with complete CH-47 produc- 
tion SAS included (with Kp 4 2) . The longitudinal control input 
C<$long-) designed in an analogous manner to the 40 knot case with 
two differences. First, the amplitude of both the low frequency and 
high frequency part was determined by trial runs such that good 
excitation resulted in u and q. Secondly, the low and high fre- 
quencies were not designed exactly at the long and short period 
modes, but slightly removed to assess whether this difference 
deteriorates identification accuracy. The time history response 
of Figure 3.5 shows excellent longitudinal velocity and pitch rate 
response, while the vertical velocity response is less than 2 ft/ 
sec. This decoupling also reflects itself in the accuracy with 
which the u and q derivatives are identified. Also, notice that 
excellent decoupling exists between the lateral and longitudinal 
motions (the decoupling again a result of the SAS) . 


The simulation data shown in Figure 3.5 was used to identify 
the u, q, and <$lonG derivatives appearing in the longitu- 
dinal equation JEq. (2.27)]. The data used in the identification 
is first contaminated with the nominal random noise, scale 
factor errors, and bias errors of Table 2.4. The identified 
derivatives are shown in Table 3.21, with those denoted by an 
asterisk held fixed to zero. The percent errors shown in 
Table 3=21 indicate the excellent accuracy attainable in the 
identified derivatives with the exception of and M^. 

is identified to 11% error in the <5 C oll ma neuver described 
next; and the nonlinear analysis. Appendix D, shows that reducing 
the data sample length also improved the prediction of Z^ and 
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Figure 3.5 Time History Response of the 6 DOF LRC Derivative Model with CH-47 
Production SAS (Kp t 2) Generated for Derivative Identification 
(Hover, <$ L0NG Input) 
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Table 3.21 


Longitudinal Derivative Identification Results 
Hover from <5 LONG Input Maneuver (6 DOF LRC Simulation, 
Random Noise, Bias, and Scale Factor Errors Included) 


PARAMETER 

TRUE 

PARAMETER 

VALUE 

IDENTIFIED 

PARAMETER 

i ERROR 

X u 

-.01857 

-.01599 

13.9 

Z u 

.02097 

* 


M u 

.04449 

.04698 


X w 

.02894 

★ 1 

ifc'H'l 


-.2598 

-.0734 

71.7 

\ 

.00935 

* 

— 

\ 

-1.3422 

-1.536 

14.5 

X «L0NG 

-1.3092 

-1.3788 

5.3 

Z «L0NG 

-.37776 

-.4248 

12.6 

^LONG 

-12.58189 

-12.2008 

3.0 

Xfi C0LL 

11.2452 

* 

— 

Z «C0LL 

-96.612 

it 

— 

^COLL 

.69449 

* 



*Held fixed to zero 












M . The time history comparisons using the identified derivative 
H 

model and the simulation data are shown in Figure 3.6. The fit is 
excellent (note that the simulation data contains random, scale 
factor, and bias errors) and, as indicated previously, the vertical 
mode (A z ) shows negligible response and the simulation response 
contains mostly noise. 

3.3.2 Collective Control Input Identification Results -- Zero 
Turbulence 


In an analogous manner to the longitudinal input maneuvers 
of the last section, the collective control input was designed 
and used with the hover 6 DOF LRC simulation. The time history 
response of the 6 DOF simulation due to the collective control 
input is shown in Figure 3.7., The vertical velocity responds 
considerably, whereas longitudinal velocity and pitch rate show 
negligible excitation. Also, the lateral time histories show 
near-zero response. 


The identified derivatives are shown in Table 3.22 resulting 
from the 5 cqlL i- n P ut maneuvers > which include the w and ^qlL 
derivatives. The derivatives previously identified from the 
(S LONG i- n P ut maneuver are shown repeated in Table 3.22 and were 
held fixed to their previously identified values (except Z w ) . 

The w and ^qlL derivatives are identified with adequate 
accuracy even with random, bias, and scale factor noise contamin- 
ation. The uncoupled vertical from longitudinal motion is, thus, 
proved to be a viable procedure for identification of the linear 
longitudinal derivatives , even though the derivatives are highly 
nonlinear with u and w as exhibited by Figure 2.5. 


Figure 3.8 shows the time history comparisons using the 
identified derivative model from the ^qqlL i n P ut maneuver and 
the simulation data of Figure 3.7. The fit is excellent and, 


as indicated previously, only the vertical motion 
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Figure 3.7 Time History Response of the 6 DOF LRC Derivative Model with CH-47 Production 
SAS (K t 2) Generated for Derivative Identification (Hover, ^cqLL * n P ut J 
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Table 3.22 


Longitudinal Derivative Identification Results at Hover 
from 6cOLL Input Maneuver (Random Noise, Bias, Scale 
Factor Errors, 6 DOF LRC Simulation) 



REF. 1 
TRUE 

®L0NG INPUT MANEUVEft 

6 C0LL INPUT ^EWE* 

PARAMETER 

PARAMETER 

VALUE 

IDENTIFIED 

PARAMETER 

% ERROR 

IDENTIFIED 

PARAMETER 

t ERROR 

X u 

-.01857 

-.01599 

13.9 

** 

— 

Z u 

.02097 


— 

— 

— 

M u 

.04449 

.04698 

5.7 

** 

— 

X w 

.02894 

* 

— 

.0295 

2.1 

Z w 

-.2598 

-.0734 

71.7 

-.2887 

11.1 

\ 

.00935 

* 

— 

.00997 

6.8 

M q 

-1.3422 

-1.536 

14.5 

** 

— 

%0NG 

-1.3092 

-1.3788 

5.3 

*★ 

— 

Z 6L0NG 

-.37776 

-.4248 

12.6 

** 

— 

m «long 

-12.58189 

-12.2008 

3.0 


— 

Xfi C0LL 

11.2452 

* 

— 

10.9680 

2.5 

Zfi COLL 

-96.612 

★ 

— 

-97.056 

.46 

^COLL 

.69450 

★ 

— 

.69016 

.58 


‘Held fixed to zero 

“Held fixed to values Identified from 8 ^q N( . maneuver 
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Figure 3.8 Longitudinal Time History Match of Identified Derivative Model 

and 6 DOF LRC Derivative Noise Contaminated Simulated Data (Hover 
'-'COLL Input, Random Noise, Bias, and Scale Factor Errors) 








excited. The 
derivative X 


excitation 
MOLL Which 


of 

is 


A is mainly due to the control 

x 1 

identified for this maneuver. 


3.3.3 Collective Input Maneuver Identification Results with 
Gust Disturbance 


This section presents the identification results with random 
gust disturbance included in the simulation. The simulation, 
control input used, and procedure are identical to those of the 
last section with the exception of the addition of the gust dis- 
turbance to the simulation. The identification program does not 
account for gust in the model structure and, as such, the results 
presented show under what conditions random gust deteriorates para- 
meter identification accuracy. The results of this section 
establish the magnitude of gust tolerable during the flight test, 
such that good identification accuracy is achieved. 

Moderate and Light Longitudinal Gust . Figure 3.9 shows the 
hover 6 DOF (Ref. 1) derivative model time history response due to 
a collective input with moderate random longitudinal gust excit- 
ation (a u = 2.13 m/sec). Comparison of this figure with Figure 3.7 
reveals the sensitivity of the response to moderate turbulence 
level. As shown in Figure 3.9, both the u and q responses are 
affected by the presence of the gust. The effect on derivative 
identification accuracy is seen by referring to Table 3.23 which 
shows that for the moderate turbulence level (a = 2.13 m/sec) 
identified parameter accuracy is degraded and significantly de- 
graded for M w and ^5 colL* Also shown in Table 3.23 is separate 
identification results where the random gust in the simulation was 
reduced to an RMS value of .91m/sec, which corresponds to light- to- 
calm air. The identified parameters shown in Table 3.23 are only 
slightly changed relative to the no- turbulence case (Table 3.22) 
and are of totally acceptable accuracy. 
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Figure 3.9 Time History Response of the 6 DOF LRC Derivative Model with CH-47 
Production SAS (Kp f 2) with Moderate Longitudinal Random Wind Gust 
Disturbance (a u = . 2.13 m/sec) 
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Figure 3.9 (Continued) 












Table 3.23 

Lcongitudinal Derivative Identification Results at Hover from 
COLL Input Maneuver with Longitudinal Wind Gust Disturbance 
(Random Noise, Bias, Scale Factor Errors, 6 DOF Simulation) 
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Held fixed to values identified from 6 ^q NG maneuver 




The time history comparisons with the simulation data of 
the identified derivative models for both the moderate and light 
turbulence levels are shown in Figure 3.10 and 3.11, respectively. 
The degraded accuracy of the moderate turbulence case is reflected 
in the poor fit of Figure 3.10 and the excessive response in q 
and 0. The time history match for the light turbulence case is 
shown to be quite good, as was reflected in the good accuracy in 
identified parameters. 

Vertical and Lateral Gust (Heavy Turbulence Level) . The 
results just presented demonstrated the effect of longitudinal 
gust on identified derivative accuracy for the collective input 
maneuver. In a similar manner, this section presents the re- 
sults where a vertical gust and then a lateral gust disturbance 
are included in the simulation. 

Table 3.24 shows the identified derivative values for these 
two cases: (1) a vertical gust with heavy turbulence, 

a w = 3.05 m/sec and (2) a lateral gust with heavy turbulence, 
o y = 3.05 m/sec. The vertical gust is shown in degrade accuracy 
in the identified stability derivatives X , Z , and M 

WWW 

unacceptably, whereas the control derivatives ^scqlL anc * ^6COLL 
retain good accuracy. The lateral gust is shown to have a neg- 
ligible effect on the identified derivative accuracy over the 
no-gust condition and thus flight testing in a cross wind would 
be an acceptable testing procedure. The time history comparison 
of the simulation model and the identified derivative model for 
both the vertical gust case and the lateral gust case is shown 
in Figures 3.12 and 3.13, respectively. The degree of time 
history match is, thus, as expected based upon the identified 
derivative accuracies of Table 3.24. 

3.4 LATERAL -DIRECTIONAL RESULTS AT THE HOVER CONDITION 

This section presents the identification results for the 
lateral-directional hover condition. The simulation model in 
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Figure 3.10 Longitudinal Time History Match of Identified Derivative Model and 

6 DOF LRC Derivative Noise Contaminated Simulated Data with Moderate 
Longitudinal Random Wind Gust Disturbance (a = 2.13 m/sec, Hover 
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Figure 3.11 Longitudinal Time History Match of Identified Derivative Model an 
6 DOF LRC Derivative Noise Contaminated Simulated Data with Light 
Longitudinal Random Wind Gust Disturbance (a u =.91 m/sec. Hover) 







Table 3.24 

Longitudinal Derivative Identification Results at Hover from a ^COLL 
Input Maneuver with Vertical and Lateral Turbulence 
(Random Noise, Bias, Scale Factor Errors, 6 DOF Simulation) 
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Held fixed to values Identified from maneuver 
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Figure 3.12 Longitudinal Time History Match of Identified Derivative Model an 
6 DOF LRC Derivative Noise Contaminated Simulated Data with Heavy 
Vertical Random Wind Gust (o w = 3.05' m/sec, Hover) 
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each case is the 6 DOF LRC derivative model with the CH-47 pro- 
duction SAS (with Kp v 2) shown in Appendix B. The procedures 
used are essentially identical to those employed for the longi- 
tudinal case. Random noise, bias, and scale factor errors are 
used and the nominal values employed are shown in Table 2.4, 
and implemented as shown in Figure 2.13. Single input excit- 
ations are employed as in the longitudinal case and designed to 
contain frequencies at the natural system modes. The procedure 
employed first identifies derivatives using a lateral input 
maneuver and then the remaining derivatives are identified using 
a rudder (lateral cyclic) input. In addition, the effect of 
random wind gust on identification accuracy is presented. 

3.4.1 Lateral and Rudder Input Maneuver Results--Zero Turbu- 
lence 


Table 3.25 presents the identified lateral parameter values 
for both the lateral control input C^lat^ maneuver and the 
rudder pedal input ( ^rud^ maneuver. The resulting accuracy for 
most of the derivatives is shown to be very good with the excep- 
tion of those derivatives contributing negligibly to the time 
response. Figure 3.14 shows the time history fit using the 
derivatives identified from the lateral control input ( 5 laT^ 
maneuver and Figure 3.15 shows the time response comparison for 
the rudder pedal input C^ruD^ maneuver. The time history fits 
are shown to be excellent for each maneuver reflecting the accur- 
acy levels shown in Table 3.25. 

As discussed previously, the lateral model employs constant 
feedback gains in the structure; thus, the primed derivatives shown 
in Table 3.25 include the effects of the roll and yaw SAS. Table 
3.26 shows the aerodynamic derivatives without the SAS as deter- 
mined from Eq. (3.9). The lateral velocity derivatives are not 
affected for the hover case since there is no sideslip feedback 
at hover. 
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Lateral Derivative Identification Results at Hover 
for Both A 6LAT and 6 rud Maneuver (Random 
Noise, Bias Scale Factors, 6 DOF Simulation) 



Held fixed to zero 

**Meld fixed to values Identified from 6^. maneuver 
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Figure 3.14 Lateral Time History Comparison of Identified Derivative Model and 
6 DOF LRC Derivative Noise Contaminated Simulated Data (Hover, 

^LAT Input, Random Noise, Bias, and Scale Factor Errors) 
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Figure 3.1!> Lateral Time History Comparison of Identified Derivative Model and 
6 DOF LRC Derivative Noise Contaminated Simulated Data (Hover, 5 rud 
I nput, Random Noise, Bias and Scale Factor Errors) 
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Table 3.26 


Lateral Aerodynamic Derivatives without the SAS 
Contribution from Equation (41) using Identified 
Values of Table 3.25 (Hover) 


PARAMETER 

Ref. 1 
TRUE 

PARAMETER 

VALUE 

MAXIMUM 

LIKELIHOOD 

ESTIMATED 

VALUE* 

% ERROR 

V P 

-.501 

CM 

O 

in 

• 

1 

0.24 

S 

-.7547 

O 

o 

00 

• 

1 

6.0 

n p 

-.0595 

-.0227 

61.8 

V 

-.0552 

-.1057 

91.4 

L r 

-.0753 

-.0822 

9.1 

N r 

-.0465 

-.00047 

98.9 


*Derived from Eq. (55): K p 3 -.064, K p = -.262 


3.4.2 Lateral Results — Heavy Turbulence 


The lateral con trol input ($ la-jO maneuver was repeated with 
heavy turbulence (a v = 3.05 m/sec) in the lateral axis. The iden- 
tified parameter values are shown in Table 3.27 and, for the most 
part, have been degraded appreciably over the zero gust case. The 
lateral input control derivatives y 6lat> l <5laT’ ancl N(S LAT are 
shown to be of fair accuracy and the stability derivative N v is 
of good accuracy. Gust, in general, is found to have the smallest 
degradation in accuracy of the control derivatives and the excel- 
lent accuracy in N, r of Table 3.27 is most likely due to the large 

w 

excitation in lateral velocity as a result of the gust, thus caus- 
ing good identif iability . The time history responses (Figure 3.16) 
show a reasonable fit between the identified derivative and simu- 
lation models. 
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Table 3.27 


Lateral Derivative Identification Results at Hover for 6^^ 

with Heavy Lateral Turbulence 
(Random Noise, Bias, Scale Factors, 6 DOF Simulation) 


PARAMETER 

REF. 1 
TRUE 

PARAMETER 

VALUE 

6 LAT INPUT 

MANEUVER 

IDENTIFIED 

PARAMETER 

% ERROR 

Y v 

-.1317 

-.2533 

92.2 

L v 

-.00623 

-.01086 

74.3 

N v 

-.00103 

-.00109 

6.5 

V P 

-4.537 

-8.79 

93.9 

L P 

-1.834 

-2.157 

17.6 

N P 

-.1653 

-.1123 

32.0 

Y r 

-.3870 

* 

— 

L r 

.5517 

★ 

— 

N r 

-2.098 

* 

— 

Y <5LAT 

1.157 

1.609 

39.1 

l <slat 

.4317 

.4769 

10.5 

N <S|_AT 

.04229 

.0315 

25.4 

Y <5rU0 

-.05506 

-- 

— 

pijn 

-.06007 

— 

— 

N(S RUD 

.1988 

■ . 

— 


★ 

Held fixed to zero 
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Figure 3.16 Lateral Time History Comparison o£ Identified Derivative Model 
and 6 DOF LRC Derivative Noise Contaminated Simulated Data with 
Heavy Lateral Turbulence (o = 3.05 m/sec, Hover) 
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3.5 IDENTIFICATION RESULTS AT THE 60 KNOT SPIRAL DESCENT 

CONDITION 

The 60 knot spiral descent flight derivatives from Ref. 1 
were used to generate a nonlinear simulation model. A nonlinear 
model was used in order to correctly represent the kinematics of 
turning flight. Control inputs were designed similar to the 40 
knot case previously discussed. The results of the 60 knot 
spiral descent flight condition are summarized below for the 
longitudinal and lateral maneuvers; the simulation model is 
shown in Appendix C. 

3.5.1 Longitudinal Identification Results 

The longitudinal inputs designed for the 40 knot condition 
were also used for 60 knots. Eigenvalues were obtained from the 
LRC derivatives at the 60 knot spiral descent flight condition 
with and without the CH-47 production SAS. This condition was 
found to he the most unstahle of the three airspeeds analyzed (due 
to longitudinal divergent root); thus, the control inputs caused 
excessive divergence in the response for both longitudinal and 
lateral input maneuvers. The open- and closed-loop eigenvalues 
are shown in Table 3.28. 

The longitudinal input maneuver used for the 40 knot condi- 
tion was used with the 60 knot simulation program. This input 
was clearly unsatisfactory due to the large degree of instability 
even with SAS-on. After several input redesign attempts, it was 
found that by setting to zero the first one-half cycle of the 
low frequency input, the response in pitch remained within accep- 
table response limits. The resulting input used for the longi 
tudinal maneuver is shown in Eq. (3.10) and is within the TAGS 
rate authority. 
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Table 3.28 


Eigenvalues for the 60 Knot Spiral Descent Flight 
Condition with and without the CH-47 Production SAS 



SAS-0FF 

CH-47 PRODUCTION 
SAS-0N 

v 2 

Long. /Vert. 
Long. Divergent 
Pitch Damping 
Dutch Roll 
Roll Damping 
Spiral 

-.112 +j .281 
.5025 
-2.607 
.08 +j .379 
-1.07 
-.0406 

-.328 +j .298 
.220 

-2.336 +j 2.637 
-.218 +j .3734 
-2.24 +j .3512 
-.00185 


>LqnG = A 1 sin 3 * 43t + A 2 sin • 525t 


cm 


A 1 = 


3.6 

0.0 


t < 1.8 
t > 1.8 


A 2 


.90 6 < t < 12 

0.0 t < 6 

0.0 t > 12 


*COLL = B 1 sin 3 - 43t + B 2 sin • 525t 


~ cm i 


B i - 


3.6 6 < t < 7 . 8 

0.0 t < 6 

« A a. V ^ rt 

U.U 


B 2 = 


.90 

0.0 


t > 6 
t < 6 


3.10) 
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The time history response for the longitudinal input of 
Eq. (3.10) for the 60 knot spiral descent condition is shown in 
Figure 3.17 using both the linear and nonlinear simulation. The 
6 long i- n P ut is shown to consist of one-half cycle of the low 
frequency input to prevent divergence of the pitch response due 
to the highly unstable pitch divergent eigenvalues. Also, the 
^vCOLL i n P ut employs the reduced value for the high frequency 
input so that saturation of the TAGS actuator is prevented. 

The longitudinal derivatives are identified with reasonable 
accuracy using either the linear or nonlinear simulation data. 
The results are shown in Table 3.29. 

The time history match of the identified derivatives and 
nonlinear simulation data is shown in Figure 3.18 for the spiral 
descent condition. The excellent match further reflects the 
accuracy in the identified derivatives. Table 3.30 presents 
the identified derivatives, standard deviations, and F-ratios 
for this .condition. 

3.5.2 Lateral Identification Results 


Lateral inputs were designe d si milar to the 40 knot condition 
and resulted in large divergence in the pitch mode 9.0 m/sec 
longitudinal velocity change]. The cross-coupling was determined 
to be very weak and the longitudinal divergence is a result of 
the pitch longitudinal divergent root. Since under actual flight 
conditions the pilot normally stabilizes this divergence, it was 
decisded that a pilot math model should be used to simulate this 
situation. Although the math model representation is only an 
approximation to pilot response, it should be adequate for the 
purposes of this study (i.e. to keep the pitch mode from diverg- 
ing) . It was assumed that a pilot model which corrected for 
pitch changes would be acceptable for stabilizing the divergent 
root . 
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Figure 3.17 Time History Comparison of the 6 DOF Nonlinear and Linear Simulation for 
the Longitudinal 60 KT Spiral Descent Condition (Hover) 
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Figure 3.17 (Continued) 
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Figure 3.17 (Continued) 
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Figure 3.17 (Concluded) 
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Figure 3.18 Longitudinal Time History Match of Identified Derivative Model and 
6 DOF Nonlinear Noise Contaminated Simulation Data (60 KT Spiral 
Descent) 
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Table 3.29 


Longitudinal Identification Results for the 
60 Knot Spiral Descent Flight Condition 


PARAMETER 

REF. 1 
TRUE 

PARAMETER 

VALUE 

0ATA FROM 6 OOF 
LINEAR SIMULATION 

OATA FROM 6 OOF 
NONLINEAR SIMULATION 

MAXIMUM LIKELIHOOD 
ESTIMATED PARAMETER 

MAXIMUM LIKELIHOOD 
ESTIMATED PARAMETER 

x u 

-.01916 

-.00873 

-.00988 

z u 

-.0740 

-.1017 

-.0536 

M u 

-.01637 

-.01814 

-.01644 

x * 

.03956 

.01310 

.0102 

Zw 

-.5395 

-.5302 

-.4479 

M w 

.06089 

.04974 

.03835 

M q 

-1.701 

-1.447 

-1.393 

^LONG 

-1.4880 

-.1656 

-.4236 

Z^LONG 

-6.1572 

.02748 

-.2508 

^5 LONG 

-15.4685 

-14.8268 

14.5512 

X «C0LL 

5.6532 

4.3132 

-.0876 

Z^COLL 

-110.196 

-113.856 

-115.38 

^COLL 

6.5118 

7.4803 

6.0866 


Table 3.30 

Longitudinal Identification Results for the 
60 kt Spiral Descent Flight Condition Using 
the Nonlinear Simulation 


PARAMETER 


REF. 1 

Tpnr 

MAXIMUM LIKELIH000 

l Kuc. 

PARAMETER 

VALUE 

IDENTIFIED 

PARAMETER 

STANDARO 

DEVIATION 

PARAMETER 

F-RATIO 

-.01916 

-.00988 

.00275 

12.9 

-.0740 

-.0536 

.00855 

39.3 

-.01637 

-.01644 

.00066 

628.3 

.03955 

.0102 

.00556 

3.36 

-.5395 

-.4479 

.01101 

1654. 

.06089 

.03835 

.00144 

697. 

*■1.701 

-1.393 

.0228 

3704. 

-1.4880 

-.4236 

.34404 

1.51 . 

-6.1572 

-.2508 


.127 

-15.4685 

-14.5512 


20463. 

5.6532 

-.0876 


.0639 

-110.196 

-115.68 

.72444 

25534. 

6.5118 

6.0866 

.09724 

3894. 
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Figure 3.19 Lateral Time History Match of Identified Derivative Model 
and 6 DOF Nonlinear Noise Contaminated Simulation Data 
(60 KT Spiral Descent) 
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Identification results for the lateral inputs are shown in 
Table 3.31, with the pilot model incorporated in the simulation to 
stabilize divergent pitch mode. The time history match to the non- 
linear simulation is shown in Figure 3.19. 


Table 3.31 

Lateral Identification Results from the Nonlinear 
Simulation Including the Pilot Model 
(60 Knot Sprial) 


PARAMETER 

REF. 1 
TRUE 

PARAMETER 

VALUE 

MAXIMUM LIKELIHOOD* 

ESTIMATED 

PARAMETER 

STAHDAR0 

DEVIATION 

F-RATIO 

r ; 

-.07208 

-.05207 

.00273 

361. 


-.23648 

-.020085 

.000446 

2014. 


.02093 

.019423 

.000584 

1105. 

y p 

-1.517 

-1.697 

.068 

620. 


-1.943 

-2.070 

.00960 

46453. 

% 

-.1797 

-.2692 

.0128 

440. 

y ; 

.1188 

— 

— 

— 

k 

.53176 

.3437 

.0191 

1934. 

K 

-1.9906 

-2.329 

.0249 

8729. 

Y <5|_AT 

13.128 

13.992 

.3523 

1567. 

L 4|_AT 

16.3307 

16.8780 

.05472 

96012. 

n «lat 

1.5906 

2.2283 

.08223 

729. 

Y< 5rUD 

-.73008 

-1.5204 

.3192 

22.7 

Li ruo 

-2.3740 

-3.7047 

.06417 

3318. 

m «ruo 

7.4370 

9.1535 . 

.08819 

10732. 


*H1tJi pilot nodal used In simulation data 


The next chapter presents the flight test plan for the CH-47 
VALT research aircraft for parameter identification based on the 
results presented in this chapter. 
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IV. FLIGHT TEST PLAN FOR CH-47 PARAMETER IDENTIFICATION 


The flight test plan presented in this chapter is based on 
the operational requirements of the VALT mission and the results 
of the previously described analytical identification study. The 
recommended flight test plan is based on a number of separate 
test points which have been selected for the task of CH-47 para- 
meter identification along the VALT trajectory. These test 
points are described, and the required pilot input are specified 
in separate sections. 

4.1 FLIGHT TEST CONDITIONS 

Table 4.1 provides an overview of the recommended flight 
test plan and each test point is described in the following 
paragraphs . 


Flight Test Point (1) : Flight test point Cl) is at a hover 

condition at 200 ft altitude. This condition is selected to 
familiarize the pilot with the recommended testing procedures. 

This point serves as a useful starting point for the remainder 
of the test by providing a "bench mark" evaluation of atmospheric 
conditions, "feel of the aircraft," and instrumentation calibra- 
tion, In addition, the CH-47 is out-of- ground effect and is 
stable with the SAS on. 

Four uncoupled inputs are applied at this point. These are: 

(1) differential collective for excitation of the longitudin- 
al mode (o L0 NG' i ’ 

(2) collective for excitation of the vertical mode C^^qlL^ * 

(3) lateral cyclic for excitation of the primary roll modes 

C6 lat ); and 
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TABLE 4.1 

Summary Table of Flight Test Conditions Recommended for 
Testing (Zero Sideslip Except Where Noted) 
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Conditions to be repeated: 13 additional maneuvers, total number of maneuvers 






















TABLE 4.1 (Concluded),- 
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(4) differential cyclic for excitation of the primary yaw 

mode ( 6 rud^ - 

Flight Test Point (2) : Following flight test point (1), the 

pilot establishes a descent rate of 2.54 m/sec at 61 m altitude. 
This is achieved by flying to 30-60 m above the altitude test 
condition, and obtaining zero forward velocity and the specified 
descent rate at approximately 23 m above the altitude test point 
at which time one of the four uncoupled inputs are applied for 
approximately 15 sec. An altitude loss of approximately 125 ft 
will be encountered during the maneuver. The trim velocities can 
be established by radio ground link to the real time radar moni- 
toring system of the VALT test facility. This procedure is re- 
peated for each of the four inputs. 

Flight Test Point (3} : The next flight test point is for 

the evaluation of in-ground effect at hover. The VALT trajectory 
point most influenced by the in-ground effect is Point D of 
Figures 2.1, 2.2 and 2.3. The VALT test point is zero forward 
airspeed at 15.2 m with a descent rate of 2.54 m/sec. Since any 
complete maneuver would then be limited in time from the safety 
of flight constraint, this point is deleted from the test plan. 
Instead, flight test point (3) is at zero descent rate at 50 ft 
altitude. This test point will isolate in-ground effect and is 
stable with the SAS-on. All four uncoupled inputs -- ^ long * ^COLL’ 
^LAT’ ^RUD'~ are testec * at this point. 

Flight Test Point (4) : The next flight test point is at 

forward speed of 20 kts. This flight test point is a VALT tra- 
jectory point and the altitude is 61 m with a descent rate of 
2,54 m/sec. This condition is stable with the SAS-on; however, 
it is also highly nonlinear. As discussed in preceding chapters, 
there is a strong effect of forward velocity on angle-of-attack 
stability (M w ) and vertical velocity on speed stability (M u ) . 
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This nonlinear M uw effect may require special inputs in order 
to provide only the linear derivatives which the objectives of 
this particular program require. 

Flight Test Point (5) : Flight test condition (5) is not a 

VALT trajectory point, but is one of three flight test points 
which are included to define CH-47 stability and control coeffi- 
cients at the forward speed of 40 kts [the other points are (6) 
and (7)]. This forward speed is predicted to be unstable with 
the SAS-on and, thus, will require pilot stabilization as well 
as additional inputs for decoupled lateral response. Test point 
(5) is a zero velocity descent rate at 40 kts with two coupled 

inputs- - ^ long" ^ COLL anc * ^LAT" ^RUD" " at an of 97.5 m. 

Flight Test Point (6) : This is a VALT trajectory point 

requiring pilot stabilization and is the same as the conditions 

of (5) except that the two coupled 6 LONG" ^COLL ancl ^LAT’^RUD 
inputs are performed at a descent rate of 2.54 m/sec. 

Flight Test Point (7) : This is not a VALT trajectory point 

but does complete the flight test data at 40 kts forward speed. 
The purpose of this point is to determine the effect of high 
descent rate and, thus, validate trends established by points 
(5) and (6) . 

Flight Test Point [8) : This flight condition (airspeed = 

60 kts, 5 00 FPM descent) is predicted to be the most unstable 
VALT trajectory point. At an altitude of 61 m, the ^LONG'^COLL 
and 6 T aT"^DTm coordinated inputs are applied. 

Flight test point (8) completes the hover and descent test 
conditions. The next two points are for the climb condition. 

The two climb conditions will provide basic information as to 
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the effect of climb rate (take-off) on stability and control 
coefficients . 


Flight Test Point (9) : This flight point is conducted at 

zero forward velocity and a climb rate of 2.54 m/sec at an alti- 
tude of 15.2 m. This is a VALT take-off tracjectory point. Four 

uncoupled inputs (^lonG’ ^GOLL* ^LAT' ^RUD^ are use< ** To achieve 
this point, the pilot may simply lift off the pad and attempt to 
establish the vertical rate as soon as possible. No safety of 
flight constraint is necessary for this powered takeoff. It is 
not necessary for the steady state vertical velocity to be estab- 
lished exactly at 15.2 m. 


Flight Test Point (10) : This flight condition is: 40 kt 

airspeed, 2.54 m/sec climb and 97.5 m altitude. This is an un- 
stable point and requires the two coordinated inputs - _<s L ONG~ ^COLL 


and 


5 LAT" 6 RUD* 


This point completes the ascent phase. 


The remaining test points are perturbations about the hover 
condition of test point (2) with commanded velocity calibration 
being established by ground-radio link. These three points are 
placed last because they require significant pilot attention. 


Flight Test Point (11a) : A sideslip is established by trans- 

lating right laterally at 10 kts. 


Flight Test Point (lib) : Flight test point (lib) is the 

same as (11a) except the sideslip is established to the left with 
a lateral velocity of 10 kts. 


Flight Test Point (11c) : Flight test point (11c) is the 

same as flight test point (2) except the CH-47C is translating 
backward with a velocity of 10 kts. 
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Flight Test Point (12) : Flight test point (12) is the same 

as flight test point (8) except a steady turn rate (one minute 
turn) is established. 

These twelve points constitute the basic flight test require 
ments . It is further recommended that certain of these test 
points be repeated at least once. These VALT trajectory points 
are 2, 4, 6, 8, and 12. 

4.2 FLIGHT TEST INPUTS 

Three basic types of inputs are required to implement the 
flight test plans. These types are as follows: 

(1) single axis inputs, 

(2) coupled axis inputs, and 

(3) pilot inputs required for stabilization or decoupling. 
The application of these inputs is now discussed. 

4.2.1, Single Axis Inputs* [(1), (2), (5), (4), (9), (11)] 

Single axis inputs are defined relative to the fuselage 
reference axis about which they are applied . For example, dif- 
ferential collective (S^q^q) is considered applied about the 
fuselage y axis, differential cyclic (<5 RUD ) about the z axis, 
and collective cyclic (<5]^) about the x axis. Collective 
(6 coll-^ a PPli e d almost along the z axis. These inputs are 
desirable because they can allow good excitation of the principal 
modes which are most affected by parameters of interest. 

The flight test single axis inputs which are recommended 
for the test plan are used only at hover and low forward speed 


*Note: Numbers in parentheses define effected flight test points 


147 



C20 kts). These inputs are shown in Figure 4.1. The input con- 
sists of two parts--a rapid doublet followed by a slower, lower 
amplitude movement. For the longitudinal input, this slower stick 
movement is a doublet and for the lateral case, the slower stick 
movement is such as to return the rotorcraft to trim. The rapid 
doublet, in all inputs, excites the shorter period modes. 

The amplitudes and frequency of these inputs are specifi- 
cally chosen to minimize nonlinear rotorcraft response affects. 

All inputs are applied with the SAS-on. 

4.2.2 Coupled Axis Inputs [(5), (6), (7), (8), (10), (12)] 

Coupled axis inputs are those which are applied about two 
or more axes during the same maneuver. Such inputs optimize the 
information in cross-coupled derivatives. 

The inputs designed for the unstable 40 kt and 60 kt points 
are shown in Figure 4.2. It should be noted that the 5 lonG 
input contains only the second h - cycle of the low frequency input. 
This was established as a requirement in the previous chapter to 
prevent the response from diverging too rapidly because of the 
unstable pitch divergence root. 

In addition to the specified inputs of Figures 4.1 and 4.2. 
there is another class of inputs. These are the inputs required 
from the pilot to stabilize or decouple the helicopter responses. 
Such inputs are themselves made in addition to the inputs speci- 
fied above. Such inputs will provide useful identification data 
and can be treated accordingly in the post- flight data reduction 
phase to supplement results from the specified inputs. 
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HORIZONTAL VELOCITY/PITCH MANEUVER \ > SIDESLIP/ROLL MANEUVER 
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Figure 4.1 Flight Test Input Design, Hover and 20 Knots (4 Individual Maneuvers) 






Control Input Design Summary 
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LONGITUDINAL MANEUVER ! LATERAL/DIRECTIONAL MANEUVER 
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Figure 4.2 Flight Test Input Design, 40 Knots and 60 Knots (2 Individual Maneuvers) 






Control Input Design Summary 
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Figure 4.2 (Concluded) 












V. CONCLUSIONS AND RECOMMENDATIONS 


S.l CONCLUSIONS 


This investigation has studied the requirements for accurate 
flight test identification of the stability and control deriva- 
tives of the CH-47 VALT Research Aircraft. From the results of 
this simulation study, it is concluded that the linear quasi- 
static aerodynamic derivatives can be identified from flight data 
with excellent accuracy for many of the important derivatives. 
Furthermore, the results of this study have isolated and expanded 
upon specific areas of parameter identification of importance to 
the general rotorcraft identification problem. In addition, the 
results have led to specific conclusions pertinent to the over- 
all VALT program system design studies. 

The principal conclusions which have led to a flight test 
specification are as follows: 

(1) Uncoupled longitudinal/lateral-directional linear 
models may be used to adequately describe the small 
perturbation responses of the CH-47 VALT helicopter 
about specific trim points. This assumes that the SAS 
is on and, therefore, longitudinal modes are decoupled 
from lateral modes. Rotor degrees of freedom have a 
minor influence on identification accuracy except for 
some high rate inputs which couple rotor/fuselage mo- 
tions due to roll rate feedback. 

(2) The representation of the landing trajectory by piece- 
wise linear models is, in general, not valid for low- 
speed perturbations from nominal. For this reason, it 
is concluded that careful attention must be made of 
possible nonlinearities at such low speeds. This may 
be performed using model structure determination pro- 
cedures and/or specifically designed inputs. 

(3) Error characteristics of the CH-47 VALT instrumentation 
and LRC data acquisition systems are acceptable for the 
task of parameter identification. 
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(4) For the cases analyzed, it is concluded that multi- 
cycle inputs offer marginal improvement in estimate 
accuracy over simple single cycle inputs such as 
doublets. Hover and 20 kt conditions require inde- 
pendent excitation of longitudinal/pitch and vertical 
modes due to nonlinearities (SAS-on permits such de- 
coupling). Higher speed inputs require simultaneous 
inputs about two axes. 

(5) Wind gusts of root-mean- square (RMS) velocity 4.14 kts in 
forward and vertical directions tend to degrade accuracy. 
Below 1.78 kts RMS velocity, negligible parameter esti- 
mate accuracy degradation is predicted. A lateral gust 
of 5.92 kts may be encountered without affecting longi- 
tudinal identification accuracy, but lateral -directional 
identification accuracy is affected. 

The results of this study have been based on a global, non- 
linear helicopter analysis (C - 81) and linear derivative models 
taken from Ref. 1. Thus, the results of the flight test program 
can be used to validate these computer simulations as well as the 
entire identification methodology upon which the flight plan is 
based. Twelve flight test points with a total of forty-nine maneu- 
vers are specified for the flight plan. Such a flight test should 
produce a most valuable data base for not only the VALT program, 
but also for rotorcraft parameter identification research and 
development in general. 

Although the emphasis in this work has been directed toward 
producing a detailed flight test specification for the CH-47, 
several issues have emerged which are of general interest to the 
overall VALT system design studies . These issues are: 

(1) The linear 10 DOF rotor/fuselage derivative model 
derived from the C-81 nonlinear simulation program 
was found to predict rotor instability with high 
gain feedback of the CH-47 production SAS. The 
present CH-47 production SAS roll rate feedback 
gain (K = -5.0 in/rad/sec ) was shown to reduce 

rotor damping by a factor of two. A roll rate 
feedback gain value of -20.0 was found to yield 
approximately 10% critical damping. (Actual 
flight tests of the CH-47 have shown that roll 
rate feedback of -17.17 produces 10% critical 
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damping.) It is concluded that the C-81 derived 
linear model is, thus, a valid and accurate tool 
for investigation of high gain feedback on rotor 
stability, and is useful for assessment of VALT 
program control system designs. 

(2) Aerodynamic nonlinearities in the landing approach 
were found to be very important in this investigation. 
Perturbations in longitudinal or vertical velocity of 
3.05 m/sec (10 ft/sec) from a nominal trim condition 
produced significant nonlinearity. As a result, linear 
theory used for navigation and control systems designs 
should only be used as preliminary to the final designs. 
More emphasis (than is currently used in rotorcraft 
control design) should be devoted to the nonlinear sim- 
ulation for early evaluation (and redesign as required) 
of linear designs and may, in fact, be necessary to 
achieve adequate vehicle performance. A time-varying 
coefficient model with coefficients varying with air- 
speed is not acceptable for modeling the CH-47 in a 
landing approach since nonlinearity also occurs with 
vertical velocity perturbations. Linear designs based 

.upon optimal techniques may, in fact, be far from opti- 
mal due to the significant nonlinear nature of heli- 
copters in a landing approach. 

(3) Random wind gust was found to produce significant response 
of the CH-47. 


5.2 RECOMMENDATIONS 


The techniques demonstrated for this effort can be extended to 
provide data for other VALT objectives. Specific items include: 

(1) Evaluation of potential rotor stability problems of 

the VALT program adaptive control system design . The 
C-81 linear rotor/fuselage model has been shown to 
accurately duplicate the feedback gain values resulting 
in actual CH-47 rotor instability trends. The present 
CH-47 production SAS roll rate feedback gain value is 
shown to reduce rotor damping by a factor of two. It 
is, thus, recommended to evaluate the VALT program 
adaptive control system designs for potential rotor 
stability problems using the C-81 rotor/fuselage linear 
simulation model. 
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( 2 ) 


Further analysis of wind shear effects on rotorcraft 
landing trajectories is a primary requirement. It is 
recognized that a logarithmic distribution is a realis 
tic assumption about the low altitude atmospheric 
description. Since low-speed rotorcraft response is 
highly dependent on this distribution, further quan- 
tification of its effect would impact near terminal 
navigation, guidance, and display functions. 

(3) Analysis of piloting techniques for parameter identi- 
fication objectives has a more general application to 
development of piloting techniques for many types of 
VTOL aircraft. In particular, the results of this 
study have shown that the CH-47 aircraft is particu- 
larly sensitive to nonlinear responses at certain 
points on the landing trajectory. Such sensitivity 
can reduce saf ety-of-flight margins. Since commercial 
VTOL operation would rely on pilots of varying experi 
ence, a training procedure, based on the types of con- 
siderations discussed in this report, could enhance 
these margins. 
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APPENDIX A 


C- 81 LINEAR DERIVATIVE MODELS (40 KTS) 


The linear models employed in this report at the 40 kt 
condition were derived from the C-81 nonlinear simulation pro- 
gram. The development of the various models, state variable 
definitions. Stability Augmentation System (SAS) , and usage of 
these models is presented in the second chapter of this report. 
The state variable representation is given here since this is 
the form used for identification, for time history simulations, 
and for computation of eigenvalues. The state variable form is 
given in Eq. (A.l) and the F and G matrices are presented 
in this appendix for each model. The following models are 
presented: 

• 10 DOF rotor/fuselage model (SAS-off) 

• First order rotor model (SAS-off) 

• 6 DOF quasi-static model (SAS-off) 

• 10 DOF rotor/fuselage model (with CH-47 production SAS) 

• 6 DOF quasi-static model (with CH-47 production SAS) 

x = Fx + Gu (A.l) 

A.l 10 DOF ROTOR/FUSELAGE MODEL (SAS-OFF) 


The state vector definition for the 10 DOF model is shown 
q. (A, 2) and the control vector is shown in Eq. (A. 3). 


T 


[9 ^ uwqvpr ^-ImR b lMR a ^TR • 


b lTR a lMR b lMR a lTR b lTR^ 


(A. 2) 
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T 

u + ^ 6 coll 6 long 5 lat 5 rud^ 


(A. 3) 


The F-matrix of stability derivatives is shown in Table A.l 
and the G-matrix of control derivatives is shown in Table A. 2. 


Table A.l 

F-Matrix (10 DOF Model, SAS-Off) 



1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

4 

- 32.06200 

0.00000 

- 0.06927 

- 0.00099 

0.06872 

- 14.58400 

0.00764 

0.00970 

- 0.17362 

5 

- 2.67600 

0.06927 

0.00000 

- 0.01191 

- 0.66346 

70.22100 

- 0.04285 

- 0.07380 

- 0.47033 

6 

0.00000 

0.00000 

0.00000 

- 0.00301 

0.00289 

- 1.15940 

0.00144 

0.14736 

- 0.18855 

7 

0.00000 

32.06200 

2.67600 

0.00103 

0.00272 

- 0.02234 

- 0.03601 

13.97000 

- 68.01300 

8 

0.00000 

0.00000 

0.00000 

0.00022 

0.00216 

- 0.03725 

- 0.00353 

0.00475 

0.00618 

9 

0.00000 

0.00000 

0.00000 

- 0.00003 

- 0.00008 

0.11857 

- 0.00041 

0.04100 

- 0.02032 

10 

0.00000 

0.00000 

0.00000 

0.04023 

0.97566 

- 56.08500 

- 0.33053 

- 50.14100 

- 12.81400 

11 

0.00000 

0.00000 

0.00000 

- 0.23282 

- 0.20557 

56.92000 

- 0.37892 

- 38.55100 

- 12.67800 

12 

0.00000 

0.00000 

0.00000 

0.13727 

0.94175 

- 19.79200 

0.11981 

51.05100 

- 0.87678 

13 

0.00000 

0.00000 

0.00000 

- 0.25110 

- 0.19747 

50.91000 

0.36920 

40.04580 

- 4.14680 

14 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

15 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0. 00000 

0. 00000 

0.00000 

16 

0 . 00000 ' 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

17 

0 . 00000 ' 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


10 

11 

12 

13 

14 

15 


16 

17 


1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

O.QOOOO 

3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

4 

0.01142 

0.17427 

0.01258 

0.20277 

- 13.12300 

- 0.82123 

- 14.97800 

- 0.62435 

5 

0.03770 

- 0.54848 

0.05659 

- 0.55953 

2.21710 

2.30340 

3.86100 

3.11020 

6 

- 0.00338 

0.05107 

0.00379 

- 0.06467 

0.98582 

- 0.21020 

2.00440 

0.32833 

7 

0.12898 

0.03271 

- 0.20399 

- 0.02745 

- 0.99561 

10.82800 

1.02250 

- 13.09200 

8 

0.02633 

0.00118 

- 0.06365 

- 0.00916 

0.26286 

6.93220 

0.36128 

- 8.03230 

9 

0.03337 

0.02084 

- 0.00641 

- 0.01704 

- 1.30750 

1.58360 

1.16900 

0.76587 

10 

- 36.61400 

- 49.52300 

- 0.00379 

0.06467 

- 39.15100 

- 937.36000 

- 2.00440 

- 0.32833 

11 

51.61800 

- 36.23900 

0.06387 

0.01172 

905.41000 

- 17.74000 

- 0.53970 

7.81370 

12 

0.00338 

- 0.05107 

- 36. 55100 

- 49.46500 

- 0.98582 

0.21020 

- 43.51500 

- 934.47000 

13 

0.02859 

0.00263 

51.51700 

- 36.02100 

0.17102 

7.02590 

898.25000 

- 15.55800 

14 

1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

. 0.00000 

0.00000 

15 

0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

16 

0.00000 

• 0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

17 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 
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Table A. 2 

G-Matrix (10 DOF Model, SAS-Off) 



1 

2 

3 

4 

1 

0.00000 

0.00000 

0.00000 

0.00000 

2 

0.00000 

0.00000 

0.00000 

0.00000 

3 

0.00000 

0.00000 

0.00000 

0.00000 

4 

1.05360 

- 0.17740 

- 0.00113 

0.07712 

5 

- 11.35200 

0.78872 

0.00258 

- 0.29688 

6 

0.02863 

- 0.34946 

0.01973 

- 0.00142 

7 

- 0.11943 

- 0.12247 

0.15600 

- 0.13777 

8 

- 0.10513 

- 0.07688 

- 0.02340 

- 0.13063 

9 

0.01821 

- 0.03622 

0.00851 

0.02279 

10 

11.20500 

4.01230 

42.83500 

71.08900 

11 

- 7.23410 

13.33200 

0.86260 

1.52850 

12 

0.64838 

0.71489 

- 31,23800 

52.84800 

13 

- 8.44060 

14.13500 

0.37912 

- 0.43187 

14 

0.00000 

0.00000 

0.00000 

0.00000 

15 

0.00000 

0.00000 

0.00000 

0.00000 

16 

0.00000 

0.00000 

0.00000 

0.00000 

17 

0. 00000 

0. 00000 

0.00000 

0.00000 


A. 2 6 DOF QUASI-STATIC MODEL (SAS-OFF) 

The state vector definition for the 6 DOF' quasi-static model 
is shown in Eq. (A. 4) and the control vector is shown in Eq . (A. 5). 
This model is derived from the 10 DOF rotor/fuselage model and the 
development is presented in the second chapter. 

T 

X = [0 <p ip u w q v p r] 

T 

U = t 6 COLL 6 LO nG 6 LAT 6 RUd1 

The F-matrix of stability derivatives is shown in Table A. 3 
and the G-matrix of control derivatives is shown in Table A. 4. 
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Table A. 3 

F-Matrix (6 DOF Quasi-Static Model, SAS-Off) 



1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

4 

- 32.06200 

0.00000 

- 0.06927 

- 0.00869 

0.06069 

- 12.84100 

0.00855 

0.13540 

- 0.41191 

5 

- 2.67600 

0.06927 

0.00000 

- 0.11703 

- 0.65657 

69.67300 

- 0.04391 

- 0.10234 

- 0.45820 

6 

0.00000 

0.00000 

0.00000 

- 0.00215 

0.00370 

- 1.32970 

0.00115 

0.13047 

- 0.16296 

7 

0.00000 

32.06200 

2.67600 

- 0.00033 

0.00082 

- 0.39209 

- 0.04274 

12.54700 

- 68.16100 

3 

0.00000 

0.00000 

0.00000 

- 0.00046 

0.00143 

- 0.32134 

- 0.00731 

- 0.83710 

- 0.07788 

9 

0.00000 

0.00000 

0.00000 

0.00010 

0.00226 

0.03148 

- 0.00189 

- 0.10754 

- 0.05625 


Table A. 4 

G-Matrix (6 DOF Quasi-Static Model, SAS-Off) 



1 

2 

3 

4 

1 

0.00000 

0.00000 

0.00000 

0.00000 

2 

0.00000 

0.00000 

0.00000 

0.00000 

3 

0.00000 

0.00000 

0.00000 

0.00000 

4 

0.79662 

0.24571 

- 0.00052 

- 0.02440 

5 

- 11.27000 

0.71150 

- 0.00061 

0.05546 

6 

0.05303 

- 0.39600 

- 0.00332 

0.00225 

7 

0.00407 

- 0.09048 

1.08590 

- 0.05503 

8 

- 0-02129 

- 0.06444 

0.56090 

- 0.05879 

9 

0.03689 

- 0.02651 

0.05345 

0.18811 


A. 3 10 DOF ROTOR/ FUSELAGE MODEL (WITH CH-47 PRODUCTION SAS) 


The state vector definition for the 10 DOF rotor/fuselage 
model with SAS is shown in Eq. (A. 6) and the control vector is 
shown in Eq. (A. 8). The development of this model is presented 


second 

chapter . 



X T = 

T ft rh i|j u 

1. *' T 

w a v t> r ai,™ bi,,*, ai^ n 

1’IIN - i*ll\ — * X IN. 



• 

b lTR a lMR 

^IMR a lTR kl.TR V 1 v 2 v 3 



V 4 V 5 V 6 

v 7 v 8 ] 

(A. 6) 

T 

u = 

[<5 coll 5 long 6 lat 5 rud^ 

(A. 7) 
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The F-matrix of stability derivatives is shown in Table A. 5 
and the G-matrix of control derivatives is shown in Table A. 6. 

The F-matrix presented in Table A. 5 includes the CH-47 pro- 
duction SAS. The recommended SAS modification (discussed in the 
second chapter) is to divide the roll rate feedback gain by a 
factor of two. This change modifies the F-matrix element occur- 
ring in the 21st row and 8th column (^21 8^ ^ rom a value of 
- 106 . 4 to -53.2. 

A. 4 6 DOF QUASI-STATIC MODEL (WITH CH-47 PRODUCTION SAS) 

The state vector definition for the 6 DOF quasi-static model 
with SAS is shown in Eq. (A. 8) and the control vector is shown in 
Eq. (A. 9). The development of this model is presented in the 
second chapter. 

T 

X = [0 <j> ip u w q v p r v^ v^ 

v 4 v 5 v 6 v 7 v 8) (A ' 8) 

“ T * ^COLL “LONG 5 LAT S RUD^ (A. 9!) 

The F-matrix of stability derivatives is shown in Table A. 7 
and the G-matrix of control derivaties is shown in Table A. 8. 

The F-matrix presented in Table A. 7 includes the CH-47 pro- 
duction SAS. The recommended SAS modification Cdiscussed in the 
second chapter) is to divide the roll rate feedback gain by a 
factor of two. This change modifies the F-matrix element occur- 
ring in the 13th row and 8th column CF^ g) from a value of -106.4 
to -53.2. 
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T ab 1 e A . 5 

F-Matrix (10 DOF Rotor/Fuselage Model with CH-47 Production SAS) 
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Table A. 5 
(Concluded) 
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Table A. 6 


G-Matrix (10 DOF Rotor/Fuselage Model 
with CH-47 Production SAS) 



1 

2 

3 

4 

1 

0.00000 

0.00000 

0.00000 

0.00000 

2 

0.00000 

0.00000 

0.00000 

0. 00000 

3 

0.00000 

0.00000 

0.00000 

0. 00000 

4 

1.05360 

- 0.17740 

- 0.00113 

0.07712 

5 

- 11.35200 

0.78872 

0.00258 

- 0.29688 

6 

0.02863 

- 0.34946 

0.01973 

- 0.00142 

7 

- 0.11943 

- 0.12247 

0.15600 

- 0.13777 

8 ' 

- 0.10513 

- 0.07688 

- 0.02340 

- 0. 13063 

9 

0.01821 

- 0.03622 

0.00851 

0.02279 

10 

11.20500 

4.01230 

42.83500 

71.08900 

11 

- 7.23410 

13.33200 

0.86260 

1.52850 

12 

0.64838 

0.71489 

- 31.23800 

52.84800 

13 

- 8.44060 

14. 13500 

0.37912 

- 0.43187 

14 

0.00000 

0.00000 

0.00000 

0.00000 

15 

0.00000 

0.00000 

0.00000 

0.00000 

16 

0.00000 

0.00000 

0.00000 

0.00000 

17 

0.00000 

0.00000 

0.00000 

0.00000 

18 

0.00000 

0.00000 

0.00000 

0.00000 

19 

0.00000 

0.00000 

0.00000 

0.00000 

20 

0.00000 

0.00000 

0.00000 

0.00000 

21 

0.00000 

0.00000 

0.00000 

0.00000 

22 

0.00000 

0.00000 

0.00000 

0.00000 

23 

0.00000 

0.00000 

0.00000 

0.00000 

24 

0.00000 

0.00000 

0.00000 

0.00000 

25 

0.00000 

0.00000 

0.00000 

0.00000 


The DOF models discussed in this report used in the simulation 
are generated from the F- abd G-matrices shown in Tables A. 7 and 
A. 8 by setting to zero all cross- coupled derivatives, thus yield- 
ing either a 3 DOF longitudinal model or 3 DOF lateral model with 
the corresponding SAS. 


A- 8 




Table A. 7 

F-Matrix (6 DOF Quasi-Static Model with CH-47 Production SAS) 
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Table A. 8 

G-Matrix (6 DOF Quasi-Static Model 
with CH-47 Production SAS) 



1 

2 

3 

4 

1 

0.00000 

0.00000 

0.00000 

0.00000 

2 

0.00000 

0.00000 

0.00000 

0.00000 

3 

0.00000 

0.00000 

0.00000 

0.00000 

4 

0.93711 

- 0.10910 

0.00012 

0.00007 

5 

- 8.05100 

- 0.03148 

- 0.00110 

0.00005 

6 

0.01764 

- 0.31958 

0.00000 

0.00021 

7 

0.05696 

- 0.01113 

1.15720 

- 0.05506 

8 

- 0.01442 

0.00538 

0.43174 

- 0.06077 

9 

- 0.00059 

- 0.04252 

0.04229 

0.19881 

10 

0.00000 

0.00000 

0.00000 

0.00000 

11 

0.00000 

0.00000 

0,00000 

0.00000 

12 

0.00000 

0.00000 

0.00000 

0.00000 

13 

0.00000 

0.00000 

0.00000 

0.00000 

14 

0.00000 

0.00000 

0.00000 

0.00000 

15 

0.00000 

0.00000 

0.00000 

0.00000 

16 

0.00000 

0.00000 

0.00000 

0.00000 

17 

0.00000 

0.00000 

0.00000 

0.00000 
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APPENDIX B 


DERIVATIVE MODELS FROM REF. 1 (HOVER) 


The linear models employed in this report at the hover 
condition were supplied by LRC. Details of the simulation 
procedure employed using these models are presented in the 
second chapter. The numerical values of the state variable rep- 
resentations for the 6 DOF quasi-static models are presented 
in this appendix for both the SAS-off and SAS-on models. The 
following two models are presented: 

• 6 DOF quasi-static model (SAS-off) 

• 6 DOF quasi-static model (with CH-47 production SAS) 

The state variable form used is shown in Eq. (B.l). 

x = Fx + Gu . (B.l) 

B.l 6 DOF QUASI -STATIC MODEL (SAS-OFF) 


The state vector definition is shown in Eq. (B.2) and the 
control vector is shown in Eq. (B.3). 


x 1 = [0 <|) ^ u w q v p r] 


u t 6 COLL ^LONG 6 LAT 6 RUD^ 


(B.2) 

(B.3) 


The F-matrix of stability derivatives is shown in 
Table B.i and the G-matrix of control derivatives is shown 
in Table B.2. 


B-l 



Table B.l 

F-Matrix (6 DOF Quasi-Static Model, SAS-Off) 



1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.11590 

3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00670 

4 

- 32.20000 

0.00000 

0.00000 

- 0.01857 

0.02894 

- 5.76930 

- 0.00127 

- 0.02869 

- 0.08620 

5 

- 3.70650 

0.00000 

0.00000 

0.02097 

- 0.25982 

0.40010 

0.00975 

0.35800 

0.43277 

6 

0.00000 

0.00000 

0.00000 

0.01356 

0.00285 

- 1.34220 

0.00027 

0.08761 

0.01388 

7 

- 3.70650 

31.98500 

0.00000 

- 0.00001 

0.00244 

0.00500 

- 0.13179 

- 10.10900 

- 0.18117 

8 

0.00000 

0.00000 

0.00000 

- 0.00005 

0.00060 

0.03561 

- 0.00623 

- 0.75473 

- 0.07538 

9 

0.00000 

0.00000 

0.00000 

0.00022 

0.00022 

- 0.14434 

- 0.00103 

- 0.05985 

- 0.04658 


Table B.2 

G-Matrix (6 DOF Quasi-Static Model, SAS-Off) 



1 

2 

3 

4 

1 

0.00000 

0.00000 

0.00000 

0.00000 

2 

0.00000 

0.00000 

0.00000 

0.00000 

3 

0.00000 

0.00000 

0.00000 

0.00000 

4 

0.93710 

- 0.10910 

0.00012 

0.00007 

5 

- 8.05100 

- 0.03148 

- 0.00110 

0.00005 

6 

0.01764 

- 0.31958 

0.00000 

0.00021 

7 

0.05696 

- 0.01113 

1.15720 

- 0.05506 

8 

- 0.01442 

0.00538 

0.43174 

- 0.06077 

9 

- 0.00059 

- 0.04252 

0.04229 

0.19881 


B.2 6 DOF QUASI-STATIC MODEL (WITH CH-47 PRODUCTION SAS) 


The state vector definition for the 6 DOF quasi-static 
model with SAS is shown in Eq. (B.4) and the control vector is 


shown in Eq. 

(B . 5) . The 

development 

of this model 

is presented 

in the second chapter. 

— 



X T 

= [6 1)7 

u w q v 

P r v x v 2 

V 3 


v 4 V S 

v 6 v ^ v 8 ] 


(B.4) 

T 

u 

= [ 6 coll 

6 LONG S LAT 

5 RUD^ 

(B.5) 
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The F-matrix of stability derivatives is shown in 
Table B.3 and the G-matrix of control derivatives is shown 
in Table B.4. 

The F-matrix presented in Table B.3 includes the CH-47 
production SAS. The recommended SAS modifications (discussed 
in the second chapter) is to divide the roll rate feedback gain 
by a factor of two. This change modifies the F-matrix element oc- 
curring in the 13th row and 8th column (F-^ 8-^ f rom a value of 
-106.4 to -53. 2. 

The 3 DOF models discussed in this report are generated from 
the F- and G-matrices shown above by setting to zero all cross- 
coupled derivatives, thus yielding a 3 DOF longitudinal model and 
a 3 DOF lateral model with the corresponding SAS. 



Table B.3 

F-Matrix (6 DOF Quasi-Static Model with CH-47 Production SAS) 
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Table B.4 


G-Matrix (6 DOF Quasi-Static Model 
with CH-47 Production SAS) 



1 

2 

3 

4 

1 

0.00000 

0.00000 

0.00000 

0.00000 

2 

0.00000 

0.00000 

0.00000 

0.00000 

3 

0.00000 

0.00000 

0.00000 

0.00000 

4 

0.79662 

0.24571 

- 0.00052 

- 0.02440 

5 

- 11.27000 

0.71150 

- 0.00061 

0.05546 

6 

0.05303 

- 0.39600 

- 0.00332 

0.00225 

7 

0.00407 

- 0.09048 

1.08590 

- 0.05503 

8 

- 0.02129 

- 0.06444 

0.56090 

- 0.05880 

9 

0.03689 

- 0.02651 

0.05345 

0.18808 

10 

0.00000 

0.00000 

0.00000 

0.00000 

11 

0.00000 

0.00000 

0.00000 

0.00000 

12 

0.00000 

0.00000 

0.00000 

0.00000 

13 

0.00000 

0.00000 

0.00000 

0.00000 

14 

0.00000 

0.00000 

0.00000 

0.00000 

15 

0.00000 

0.00000 

0.00000 

0.00000 

16 

0.00000 

0.00000 

0.00000 

0.00000 

17 

0.00000 

0.00000 

0.00000 

0.00000 
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APPENDIX C 


NONLINEAR SIMULATION MODEL 


The nonlinear simulation model was derived by representing 
the fixed operating point linear derivative model by nonlinear 
polynomials. The polynomials were developed from tabular 
derivative data as a function of airspeed and descent rate, 
obtained from Ref. 1 (for the hover condition and from LRC 
(for the 60 kt spiral descent condition) . The structure of the 
polynomials used in the curve fit are 

f i * k l + k 2 x 4 * k 3 x 4 * k 4 x 5 * k 5 x 5 + k 6 x 4 x 5 

♦ kjX^Xj +k 8 x^x 5 +k 9 xjxj (C.l) 

where x 4 = u and x 5 = w. The nonlinear simulation model 
(with SAS-off) is shown in Eq. (C.2). 

Tables C.l and C.2 present the numerical values of the 
polynomial coefficients for the hover and 60 kt spiral descent 
conditions, respectively. 


C-l 




C-2 



Polynomial Coefficients Used for the Hover Nonlinear Simulation Model 
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APPENDIX D 


NONLINEAR DERIVATIVE MODEL 


Input design studies were based on linear derivative models 
for simulation of the CH-47. The evaluation of derivative non- 
linearity was limited to a qualitative assessment of the im- 
portance of nonlinearity on derivative identification accuracy. 

This qualitative assessment was based upon examination of the 
primary speed dependent derivatives and the nature of their varia- 
tions about the trim condition. The specific conclusions reached, 
based upon this qualitative assessment, are that: (1) for airspeeds 

of 20 kts and below, longitudinal maneuvers should be designed to 
independently excite airspeed while holding descent rate near its 
trim value, and the reverse procedure (constant airspeed, variable 
descent rate); and (2) for airspeed of 40 kts and above, nonlinearity 
is not significant and both u and w can be excited during one man- 
euver. The present discussion is aimed at validating this as- 
sumption and quantifying the resulting derivative identification 
accuracy using a nonlinear simulation model. 

The nonlinear simulation model used throughout this analy- 
sis is derived from linear derivatives obtained at various trim 
speeds and descent rates from the LRC CH-47 nonlinear simulation 
program HELICOP Ref. [1]. A polynomial is selected to represent 
the variation in the derivative values as a function of u and w 
and a least square curve fit is used to determine the coeffi- 
cients of the polynomial. Equation (D.l) shows the polynomial 
form selected for representing the derivative variation with u 
and w. This is a second-order polynomial with all possible 
cross terms between u and w. This model was selected based 
upon examination of the variation in the derivatives with u 
and w and is valid for airspeed changes of _+20 kts and ver- 
tical speed changes of +^500 FPM. 



(D.l) 




2 2 

k^ + k2U + kjU + k^w + k^w + k^uw 

2 2 2 2 

+ kyuw + kgU w + kgU w 


where i = 1,2,3, ... ,60 represents a different set of coeffi- 
cients k.| for each stability derivative, and u and w are 
deviations from a reference trim (values for K . ' s may be found in 
Appendix C) . 


As an example, Figure 2.5 shows the derivatives M u and M w 
as a function of airspeed and descent rate. The polynomial of 
Eq. (D.l) for M u at the hover condition is derived by selecting 
nine numerical values of the derivative M u corresponding to the 
nine values of (u,w) coordinates (see Figure 2.5). This estab- 
lishes nine equations of the form of Eq. (D.l) in nine unknowns 
(the K's). The equations are solved for the unknown K’s and 
the resulting polynomial fit passes uniquely through the nine 
points used in the fit. 

The structure of the nonlinear simulation used throughout 
this analysis is shown in Eq. (D.2). 


” • 

X 


F (u,w) 

G ( u , w ) 


X 

+ 

G (u,w) 

• 

V 


B c K c 

w 

A c 


V 


0 


where the x vector represents the rigid body state and the 
v vector represents the SAS state variables previously defined. 
The F and G matrices are nonlinear polynomial functions of 
the state variables u and w. An alternate way of viewing 
Eq. (D.2) is to consider the F and G matrices as functions 
of time, since both u and w are, in fact, time functions. 

The computer simulation allows monitoring the time variation of 
F and G along with the time history response of the state 
variables . - 
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The full nonlinear model used at LRC from which the poly- 
nomial representation was derived is more accurate since the 
external force-moment components are nonlinear with the other 
state variables v, q, p, r, and control inputs. However, 
over the +_20 kt airspeed about a reference trim, the significant 
nonlinearity is due to u and w and thus Eq.(D.2) provides 
an accurate representation of the Ref. 1 simulation data. 

To further the understanding of the nonlinear polynomial 
simulation, a time history simulation is presented and discussed. 
It was previously found that changes in w must be kept small 
in order to identify the linear derivative M^; otherwise, non- 
linearity would be significant. Thus, at hover a longitudinal 
input 5 long WaS usec * t0 exc ^- te u > anc * 9 states while 

w remained near its trim value. This maneuver is repeated here 
using the nonlinear simulation and is shown in Figure D.l. The 
linear simulation is also shown by the solid line superimposed 
upon the nonlinear response. The nonlinear time response of u, 
q, and 0 shows excellent comparion to that of the linear 
model except near the end of the data record where t > 13 sec. 
Also, it is seen that w shows much larger discrepancy between 
linear and nonlinear response indicating the possibility of 
nonlinearity in the Z-force equation. Although the response 
comparison of Figure D. 1 gives a quantitative measure of the 
degree of nonlinearity, it provides little information as to 
the actual variation in the derivatives. 

The variation with time of the stability and control deriva- 
tives are an output from the nonlinear C-81 simulation. Figure 
D.2 shows this variation for the longitudinal stability and 
control derivatives for this maneuver. This figure permits 
immediate interpretation of the nonlinear behavior of the sys- 
tem. It is emphasized that at this condition (hover), the 5 lqNG 
input was chosen to retain linearity in the important derivative 
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Figure; D.l Time History Comparison of the 6 DOF Nonlinear and Linear 
Simulation for the 6 inwr Input Maneuver (Hover, SAS-On) 
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Figure D.l (Concluded) 








Figure 1).2 Time Variation of Longitudinal Derivatives for the 
Maneuver (Nonlinear Simulation, Hover) 















M u - is highly nonlinear with u and its effect in the • 

response is not significant since w is maintained near trim. 
The derivatives M u , and are near ^y constant and 

equal to their trim value. The derivatives X u> and 

^ 5 LONG are i- cant ly nonlinear. The source of nonlinearity 

is mainly attributed to large excursions in u. 

The problem of trying to identify stability derivatives 
for a linear model when the derivatives are a function of time 
can sometimes be approached by reducing the data sample period. 
To illustrate this point, a second derivative identification 
run was made in which the data length was reduced to 12 sec 
thereby keeping airspeed constant. The original data length 
used was 15 sec and a reduction to 12 sec resulted in a 10 ft/ 
sec maximum value for u. This improved identification accuracy 
considerably. The results are shown in Table D.l. 


Table D.l 

The Effect of Nonlinearity on Identification Results 
Using the Nonlinear Simulation for the 6 
Input Maneuver (Hover) 



TRUE 

PARAMETER 

VALUE 

IDENTIFIED PARAMETER VALUE USING: 

PARAMETER 

LINEAR 
SIMULATION 
(15 SEC) 

NONLINEAR 
SIMULATION 
(15 SEC) 

NONLINEAR 
SIMULATION 
(12 SEC)* 

X u 

-.01857 

-.0158 

.00448 

-.01905 

"u 

.04449 

.04968 

.004619 

.04944 


-.2598 

-.0692 

-.7056 

-.2046 

M q 

-1.3422 

-1.522 

-1.424 

-1.476 

X* 

-1.3092 

-1.1160 

-.5988 

-1.11588 

^LONG 

-.37776 

-.5352 

-.5976 

-.54036 

^LONG 

-12.58189 

-12.1535 

-12.9606 

-12.3071 


•Reduction of data length to 12 sec removed major nonlinearity due 
to u. 
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Figure D.3 shows the time history match of the identified 
derivative model and the 6 DOF nonlinear simulation using only 
12 sec of data in the identification. The simulation data con- 
tains random noise, bias, and scale factor errors as previously 
employed in the last chapter. The reduction of data length to 
12 sec completely eliminated degradation in identified deriva- 
tives due to nonlinearity as shown in Table D.l and Figure D.3. 

The results just presented demonstrate that accurate linear 
derivatives can be identified from a nonlinear simulation with 
noise contamination in the form of random errors, bias, and 
scale factor errors representative of flight test data. 
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Figure D.3 Longitudinal Time History Match of Identified Derivative Model 

and 6 DOF Nonlinear Noise Contaminated Simulation Data Using 12 Sec 
of Data (Hover, (S^q^ Input, Random Noise, Bias, and Scale Factor 
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